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NOMENCLAfURE 


A.  ,B. ,U0^  - constants  in  ith  material  property  equation  [Eq.  (5)] 


= Young's  moduli  in  principal  material  directions 


= Young's  modulus  at  45°  to  r and  z directions 


secant  modulus 


= shear  modulus  in  rz  plane 


= constant  in  energy  weighting  function  [Eq.  (29)] 


= radial,  axial,  and  circumferential  directions  [Fig.  2] 
= strain  energy  function  [Eq.  (4)] 


= rectangular  coordinates  aligned  with  principal  material 


di recti ons 


= rectangular  coordinates  in  nonprincipal  material 


di recti ons 


= angle  of  orientation  of  off-axis  tests 


= shear  strain  in  rz  plane 


= strains  in  principal  material  directions 


= strains  in  x',  y',  z’  directions 


vy ,,  \>r^  vZQ  = Poisson's  ratios  for  principal  material  directions 


[Eq.  (2)] 


3 y G y G y T 

r z 9 rz 


= axi symmetric  stresses  in  principal  material  coordinates 


Subscripts 


c = compression 


t = tension 


w = weighted 


Superscripts 


mn  = principal  material  coordinates 


pq  = principal  stress  coordinates 
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1.  INTRODUCTION  ! 

Artificial  graphites  have  been  used  for  the  past  decade  in  reentry  vehi- 
cle nosetips.  Accurate  stress  analysis  procedures  are  essential  in  the  design 
of  nosetips  because  of  exacting  performance  requirements.  A successful  nose- 
tip  design  must  have,  on  one  hand,  enough  material  to  avoid  a burn-throunh 
failure  due  to  ablation  from  aerodynamic  heatinc.  On  the  other  hand,  a nosetip 
cannot  have  so  much  material  that  a thermal  stress -generated  failure  results 
due  to  high  thermal  gradients.  The  key  element  in  any  stress  analysis  is  the 
stress-strain  (constitutive)  relationship  or  material  model. 

The  stress  analysis  problems  inherent  to  reentry  vehicle  nosetip  design 
we  re  discussed  by  Jones  [1]  in  1967  along  with  numerous  specific  problems  by 
other  authors  in  the  same  conference  proceedings  volume.  Since  then,  nosetip 
stress  analysis  technology  has  been  periodically  reviewed,  In  1970,  results  of 
the  aborted  Graphite  Advanced  Development  Plan  were  summarized  from  a stress 
analysis  and  material  modeling  point  of  view  [2].  Also,  an  extremely  valuable 
interagency  conference  was  held  by  the  Atomic  Energy  Commission  [3].  More 
recently,  Schneider  et  al  [4]  described  much  of  the  current  thermodynamic,  and  j 

elastic  stress  analysis  procedures  for  reentry  vehicle  nosetips,  Also,  Jones  j 

and  Koenig  [5]  discussed  the  material  modeling  characteristics  necessary  for 
stress  analysis  of  reentry  vehicle  nosetips. 

j 

Graphites  used  for  reentry  vehicle  nosetips  are  manufactured  in  cylindri-  j 

cal  billets.  In  the  billet  comoaction  process,  the  more  or  less  elliptically  ! 

I 

shaped  highly  anisotropic  graphite  particles  tend  to  lie  with  their  lonq  axes  1 

! 

perpendicular  to  the  billet  axis.  Moreover,  the  particles  are  randomly  oriented  4 

in  planes  perpendicular  to  the  billet  axis.  The  properties  in  the  plane  per- 
pendicular to  the  billet  axis  are  therefore  isotropic.  However,  perpendicular 
to  that  plane,  i.e.,  parallel  to  the  billet  axis,  the  properties  are  quite 

1 


different.  Accordingly,  ATJ-S  graphite  is  a transverse!'/  isotropic  material 
with  its  principal  material  directions  aligned  with  the  cylindrical  coordinate 
system,  r,  z,  0,  of  the  billet  shown  in  Fig.  i.  Reentry  vehicle  nosetips  are 
made  as  axisymmetric  bodies  with  the  symmetry  axis  being  the  billet  axis.  Thus, 
the  billet  cylindrical  coordinate  system  is  also  used  for  a nosetip. 

Two  common  nosetip  geometries,  shell  tips  and  plug  tips,  are  shown  in  Fig. 
2.  The  internal  contour  of  the  shell  tip  must  be  designed  to  have  a wall  thick- 
ness sufficient  to  satisfy  the  ablation  requirements,  but  thin  enough  tc  avoid  a 
thermal  stress  failure.  For  the  plug  tip,  the  exposed  length  of  the  plug,  the 
shank  diameter,  and  the  interface  angle  between  the  plug  and  the  surrounding 
heat  shield  material  can  be  varied  to  obtain  a plug  tip  that  can  survive  a tra- 
jectory from  a thermal  stress  versus  ablation  point  of  view.  For  both  nosetip 
geometries,  the  thermal  and  mechanical  stresses  and  strains  cannot  be  predicted 
to  the  accuracy  required  to  obtain  the  desired  design  confidence  by  using 
previously  available  material  models. 


1.1  MATERIAL  NONLINEARITIES 

A new  phenomenon  known  as  "biaxial  (multiaxial)  softening"  was  found  in 
tests  of  ATJ-S  graphite  [5-9].  Biaxial  softening  is  characterized  by  the  dr-vel' 
opment  of  slightly  larger  strains  in  biaxial  tension  than  in  uniaxial  tension, 
as  shown  in  Fig.  3,  in  contradiction  to  what  miaht  be  anticipated  on  the  basis 
of  conventional  Poisson  effects.  This  phenomenon  is  attributed  to  plastic 
volume  changes  resulting  from  internal  tearing  or  microcracking.  Other  distin- 
guishing characteristics  of  graphites,  which  render  them  more  difficult  to 
theoretically  model  than  most  metals,  are  transversely  isotropic  nonlinear 
stress-strain  behavior;  different  tensile  and  compressive  stiffnesses  and 
strengths;  and  unloading  behavior  which  is  dependent  on  stress  level.  Thus, 
the  usual  concepts  of  classical  plasticity  are  not  applicable  to  graphites. 
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1.2  DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION 

A significant  characteristic  of  graphites  is  their  different  behavior 
under  tensile  and  compressive  loads.  Both  the  elastic  moduli  (stiffnesses) 
and  the  strengths  in  principal  material  property  directions  of  these  orthotro- 
pic materials  are  differe.-t  for  tensile  loading  than  for  compressive  loading. 
This  characteristic  behavior  is  shown  schematically  in  the  stress-strain  curve 
of  Figure  4.  This  phenomenon  is  but  one  of  several  differences  that  make  com- 
posite materials  more  difficult  to  analyze  (and  hence  design)  than  the  more 
common  structural  materials  such  as  aluminum. 

Both  fiber-reinforced  and  granular  composite  materials  have  different 
moduli  in  tension  and  compression  as  displayed  in  Table  1.  Unidirectional 
glass  fibers  in  an  epoxy  matrix  have  compressive  moduli  20%  lower  than  the 
tensile  moduli  [10],  For  some  unidirectional  boron/epoxy  fiber-reinforced 
laminae,  the  compressive  moduli  are  about  15-20%  larger  than  the  tensile  moduli 
[11].  In  contrast,  some  unidirectional  graphite/epoxy  fiber -reinforced  laminae 
have  tensile  moduli  up  to  40%  greater  than  the  compressive  moduli  [11].  Other 
fiber-reinforced  composites  such  as  carbon/carbon  have  tensile  moduli  from  two 
to  five  times  the  compressive  moduli  [12],  Thus,  no  clear  patter  of  larger 
tensile  than  compressive  moduli  or  vice  versa  exists  for  fiber-reinforced  com- 
posite materials.  A plausible  physical  explanation  for  this  puzzling  circum- 
stance has  yet  to  be  made. 

For  granular  composite  materials,  the  picture  is  no  clearer.  ZTA  graph- 
ite has  tensile  moduli  as  much  as  20%  lower  than  the  compressive  moduli  [13]. 

On  the  other  hand,  ATJ-S  graphite  has  tensile  moduli  as  much  as  20%  more  than 
the  compressive  moduli  [14]. 

Many  other  materials  have  different  tensile  and  compressive  moduli.  Which 
modulus  is  higher  may  depend  on  the  fiber  or  granule  stiffness  relative  to  the 


6 


FIGURE  4 STRESS-STRAIN  CURVE  FOR  A MATERIAL  WITH 

DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION 
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TABLE  1 

AND  COMPRESSIVE  MODULI  RELATIONSHIPS 


FIBEROUS 

REPRESENTATIVE 

OR 

MODULI 

GRANULAR 

RELATIONSHIP 

FIBEROUS 

Et  = 1 -25Ec 

FIBEROUS 

Ec  = 1 -2Et 

FIBEROUS 

Et=  l-4Ec 

FIBEROUS 

E.  = 2-5E 
t c 

GRANULAR 

E = 1.25E 
c t 

GRANULAR 

Et=1.2Ec 

GRANULAR 


LA  * 


matrix  stiffness.  Such  a relationship  would  influence  whether  the  fibers  or 
granules  tend  to  contact  and  hence  stiffen  the  composite.  A general  physical 
explanation  of  the  reasons  for  different  behavior  in  tension  and  compression 
is  not  yet  available.  Investigation  of  the  micromechanical  behavioral  aspects 
of  composite  materials  may  lead  i rational  explanation  of  this  phenomenon. 
Until  such  an  explanation  is  available,  the  apni^-nt  behavior  can  be  used  in 
analyzing  the  stress-strain  behavior  of  materials.  That  is,  even  without  know- 
ing why  the  materials  behave  as  they  do,  we  can  model  their  apparent  behavior. 

Actual  stress-strain  behavior  is  probably  not  as  simple  as  shown  in  Fig.  4 
Instead,  a nonlinear  transition  region  mav  exist  between  the  tensile  and  com- 
pressive linear  portions  of  the  stress-strain  curve  [15].  The  measurement  of 
strains  near  zero  stress  is  difficult  to  perform  accurately,  but  the  stress- 
strain  behavior  might  be  as  shown  in  Figure  5 wherein  replacement  of  the  actual 
behavior  by  a bilinear  model  is  offered  as  a simplification  of  obviously 
nonlinear  behavior.  For  most  materials,  the  mechanical  property  data  is  not 
sufficient  to  justify  use  of  a more  complex  material  model.  P possible  disad- 
vantage of  bilinear  stress-strain  curve  approximation  is  the  discontinuity  in 
slope  (modulus)  at  the  origin  of  the  stress-strain  curve. 

Given  that  the  uniaxial  stress-strain  behavior  is  approximated  by  a 
bil  near  representation,  the  definition  remains  of  the  actual  multiaxial  stress 
strain,  or  constitutive,  relations  that  are  reauired  in  structural  analysis. 
Over  the  past  ten  years,  Ambartsumyan  and  his  co-workers,  [References  16  to  1?] 
in  the  process  of  obtaining  solutions  for  stress  analysis  of  shells  and  bodies 
of  revolution,  defined  a set  of  stress-strain  relations  that  will  be  referred 
to  herein  as  the  Ambartsumyan  material  model.  Jones  [20]  applied  the  model  to 
the  problem  of  buckling  under  biaxial  loading  of  circular  cylindrical  shells 
made  of  an  isotropic  material.  However,  in  application  of  the  Ambartsumyan 
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material  model  to  orthotropic  materials,  certain  deficiencies,  such  as  a non- 
symmetric  compliance  matrix  in  the  stress-strain  relations  * , are  apparent. 

Jones  [22]  also  applied  modified  bilinear  stress-s{raii.  relations  to  the 
buckling  of  shells  with  multiple  layers  of  orthotropic  matenV*  wil.i  different 
moduli  in  tension  and  compression.  His  modifications  consisted  of  weinhtino 
tensile  and  compressive  compliances  according  to  the  proportions  of  the  princi- 
pal stresses  in  order  to  obtain  a symmetric  compliance  matrix.  Isabekian  and 
Khachat-yan  [23]  made  the  Ambartsumyan  material  model  compliance  matrix  sym- 
metric just  by  enforcing  certain  relations  between  the  material  properties. 
Their  relations  are  part  of  the  basis  for  the  present  report. 

1 .3  SUMMARY  REMARKS 

In  this  report,  a new  material  model  is  presented  for  the  initial  loading 
behavior  of  graphites.  This  inherently  simple  material  model,  which  is  based 
on  a new  deformation  theory  of  urthotropic  plasticity,  has  excellent  potential 
for  the  description  of  the  biaxial  softening  phenomenon.  A pragmatic  engineer- 
ing approach  is  used  to  analyze  an  important  problem  that  has  thus  far  not  been 
amenable  to  solution  by  classical  plasticity  approaches.  The  model  is  partic- 
ularly applicable  to  materials  for  which  a finite  plastic  volume  change  is 
exhibited;  i.e.4  the  usual  restriction  of  plasticity  theories  to  zero  plastic 
volume  change  is  discarded.  The  model  is  examined  for  sensitivity  to  (1)  a 
nonsymmetric  compliance  matrix,  (2)  how  the  material  model  constants  are  obtain- 
ed, and  (3)  variations  in  material  properties.  Finally,  favorable  comparisons 
are  made  between  predicted  strain  values  obtained  with  the  new  material  model 
incorporated  in  the  SAAS  III  finite  element  computer  program  [24]  and  experi- 
mental strain  values  from  biaxial  tests  conducted  by  Jortner  [6-9], 

This  report  is  divided  into  several  parts:  nonlinear  deformation  model, 
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different  moduli  in  tension  and  compression  elastic  model,  nonlinear  deforma- 


tion model  for  materials  with  different  moduli  in  tension  and  compression, 
comparison  of  predicted  and  measured  off-axis  uniaxial  response,  and  comparison 
of  oreuicted  and  measured  biaxial  strain  response.  This  approach  is  taken 
to  build  up  the  final  comprehensive  model  by  examining  the  individual  facets 
of  the  parts  of  the  model  separately.  We  therefore  are  able  to  assimilate 
the  basic  concepts  of  the  final  modal  in  a palatable  piecemeal  fashion. 
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2.  NONLINEAR  DEFORMAT I-'  N MODEL 


Current  material  models  are  not  adequate  for  prediction  of  multiaxial 
nonlinear  stress-strain  reponse  of  graphitic  materials.  Thus,  a new  material 
mocel  will  be  formulated  in  the  following  sections  The  general  characteris- 
tics of  the  new  model  are  described  ;r:  Section  2.1.  The  effect  of  the  con- 
:lar.is  in  the  new  material  model  is  determined  in  Section  2.2.  iteration 
and  convergence  of  the  material  model  to  actual  nonlinear  stress-strain  be- 
havior is  demonstrated  in  Section  2.3.  Finally,  the  sensitivity  of  the  model 
to  nonsynmetry  of  the  compliance  matrix,  the  values  of  the  material  constants; 
and  inherent  material  property  data  scatter  is  discussed  in  Section  2.4. 

First,  however,  the  graphite  material  model  due  to  Weng  [25,26]  will  be 
described.  He  uses  a deformation  theory  of  plasticity  approach  to  develop  a 
generalized  nonlinear  stress-strain  relation  for  a class  of  bodies  made  of 
transversely  isotropic  materials.  Two  important  assumptions  are  made  in  his 
formulation  in  addition  to  the  restriction  to  monotonically  increasing  propor- 
tional loads.  The  first  assumption  is  that  the  deformations  do  not  alter  the 
degree  or  direction  of  the  anisotropy  of  the  material.  The  implications  and 
validity  of  this  assumption  for  artificial  graphites  will  be  explored  further 
in  Section  2.4.  The  second  assumption  is  that  the  shape  of  the  nonlinear 
stress-strain  curves  is  identical  for  all  stress  components.  Consequently, 
the  Poisson's  ratios  in  Weng's  model  are  constant  for  all  deformation  levels. 
Weng's  generalized  stress-strain  relation  for  an  anisotropic  material  is 
6 

e,  = X S.  -o  -4>  i = 1,  2 6 (1) 

1 j=l  1J  J 

where  nonlinear  effects  are  introducted  through  the  variable  proportionality 
function,  $.  In  addition,  = strain  components,  S^j  = compliances,  and  0j  = 
stress  components.  The  function  9 is  the  effective  strain  divided  by  the  effec- 
tive stress  and  is  obtained  by  curve-fitting  both  uniaxial  and  biaxial  experi- 
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mental  data.  However,  for  nonisotropic  materials,  no  physically  meaningful 
definition  exists  for  effective  stress  and  effectiv  • strain. 

For  artificial  grpahites,  the  normal  stress  - normal  strain  curves  are 
similar,  but  not  identical,  in  shape  (see  Fig.  6).  Also,  the  Poisson's  ratios 
are  not  always  constant;  they  can  decrease  substantialy  due  to  microcracking 
related  to  the  biaxial  softening  phenomenon.  Accordingly,  Weng's  model  cannot 
be  used  to  obtain  engineering  correlation  between  theory  and  experiment. 


2.1  GENERAL  DESCRIPTION  OF  MATERIAL  MODEL  CHARACTERISTICS 

In  the  new  material  model,  Weng's  assumption  of  identically  shaped  stress- 
strain  curves  is  discarded  as  is  his  restriction  to  constant  Poisson's  ratios. 

One  objective  of  this  research  is  to  incorporate  the  new  material  model  in  the 
finite  element  stress  analysis  computer  program  SAAS  III  [24].  For  that  program, 
the  most  general  case  that  can  be  analyzed  is  the  axisymmetric  deformation  of 
an  axisymmetric  body  made  of  a nonlinear  orthotropic  material  and  loaded  in  an 
axisymmetric  manner.  Consequently,  the  new  material  model  is  discussed  for 
this  case.  Of  course,  the  concept  of  the  model  is  applicable  to  more  general 
loading  and  deformation,  as  well  as  more  general  materials. 

The  orthotropic  stress-strain  relations  for  an  axisymmetric  elastic  body 
under  axisymmetric  load  are: 


(2) 
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Directions  denoted  by  the  subscripts  r,  z,  and  3 are  principal  material  direc- 
tions. The  elastic  material  constants  in  the  compliance  matrix  of  Eq.  (2)  are 
Ef  = Young's  modulus  in  the  r direction 

Ez  = Young's  modulus  in  the  z direction 

Ee  = Young's  modulus  in  the  9 direction 

vrz  = -cz/er  for  the  loading  ar=o  (all  other  stresses  zero) 

vrg  = -ee/er  for  the  loading  ar=o  (all  other  stresses  zero) 

vz6  = -cV£:r  ^or  ^oac^in9  °z=a  (all  other  stresses  zero) 

Grz  = Shear  modulus  in  the  rz  plane 
The  reciprocal  relations  of  orthotropic  elasticity 

vrz/^r  = vzr/^z>  vr8^r  = ver/^©»  vze^z  = vqz^  (3) 

can  be  used  to  express  alternative  definitions  for  the  Poisson's  ratios  in 

terms  of  the  seven  independent  elastic  constants. 

In  the  present  material  model,  the  seven  independent  material  properties 
are  allowed  to  vary  (independently)  with  the  state  of  stress  and  strain.  (In 
constrast,  Weng  introduced  nonlinear  effects  by  using  a single  proportionality 
function,  <{>.)  In  order  to  relate  the  material  properties  to  a multiaxial 
stress  state,  the  energy  function, 

U = (orer  + czez  + ooe0  + Trzyrz)/2  (4) 

is  used.  This  function  is  half  the  function  used  by  Weng  in  his  material  model 
for  graphites.  Since  unloading  is  not  considered,  the  distinction  between  in- 
elasticity and  nonlinear  elasticity  is  not  relevant  at  this  point.  Consequent- 
ly, the  right-hand  side  of  Eq.  (4)  corresponds  to  the  strain  energy  of  the 
equivalent  linear  elastic  system  at  each  stage  of  deformation  in  a deformation 
theory  of  plasticity.  The  composition  of  U fo**  a uniaxial  loading  case  (ox=o 
and  all  other  stresses  zero)  is  depicted  in  Fig.  7.  Two  important  points 
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should  be  noted  about  the  function  U;  (1)  it  is  a scalar  and  (2)  its  value 
depends  only  on  the  final  state  of  stress  and  strain,  not  on  the  manner  by 
which  these  states  were  reached  [27].  Thus,  in  a computer  program,  the  strain 
energy  U is  found  by  the  multiplication  of  the  eight  values  in  Eq.  (4)  rather 
than  numerical  integration  of  the  stress-strain  histories  representing  the 
strain  energy  density  or  the  complementary  energy  density. 

In  the  present  approach,  equations  are  developed  to  relate  the  material 
properties  to  the  strain  energy  U.  As  long  as  the  material  property  changes 
monotonically  (increasing  or  decreasing)  for  increasing  values  of  the  strain 
energy,  a unique  relationship  for  the  properties  in  terms  of  U does  exist. 
Material  properties-versus-U  curves  can  be  obtained  from  uniaxial  stress-strain 
data  for  the  material  and  are  approximated  by 


Material  Property.;  = 1 - B1-(U/U0i)C'i 


where  i is  the  numerical  designation  of  the  specific  material  property.  Equa- 
tion (5)  can  be  fit  to  all  material  property-versus-U  curves  quite  effectively 
for  ATJ-S  graphite  and  should  be  flexible  enough  to  represent  many  other 
engineering  materials.  Although  more  general,  Eq.  (5)  is  similar  in  form  to 
the  well-known  Ramberg-Osgood  stress-strain  relation  [28]. 

The  constant  A in  Eq.  (5)  is  the  initial  (elastic)  value  for  the  property. 
Thus,  the  dimensions  of  the  material  property  are  reflected  in  the  dimensions 
of  A.  The  constant  U0  is  included  in  Eq.  (5)  to  establish  a dimensionless 
ratio  and,  thus,  a dimensionally  correct  equation  independent  of  the  value  of 
C.  Both  U0  and  U have  the  dimensions  of  energy  density  (stress).  The  numeri- 
cal value  of  U0  is  arbitrary  and  is  selected  before  the  values  of  B and  C are 
determined. 

Constants  B and  C in  Eq.  (5)  are  dimensionless;  their  numerical  values  can 
be  determined  by  at  least  two  different  approaches.  First  a least  squares  fit 
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(regression  analysis)  approach  can  be  used.  A second  simpler  technique  is  to 
solve  simultaneously  for  the  two  constants  using  data  at  only  three  points  on 
the  property-versus-U  curve.  Since  the  value  of  U can  be  larger  under  multi- 
axial  loading  than  under  uniaxial  loading,  the  material  property  equations  are 
often  extrapolated  for  multiaxial  cases.  For  this  reason,  the  second  approach 
is  recommended.  The  three  points  can  easily  be  selected  '.o  prejudice  the 
material  property  equacion  toward  a better  representation  of  the  material 
property-versus-U  curve  at  large  values  of  U.  Thus,  difficulties  that  would 
otherwise  arise  under  extrapolation  are  avoided. 

Typical  material  property  curves  for  ATJ-S  graphite  are  shown  in  Figure 
6.  Schematics  of  those  curves  and  associated  material  property  equations 
for  a modulus  and  a Poisson's  ratio  are  shown  in  Figures  8 and  9,  respectively. 
The  type  of  uniaxial  stress-strain  data  required  to  define  the  material  property 
equations  is  also  shown  in  these  figures.  Secant  moduli  are  used  when  the 
stress-strain  curve  ceases  to  be  linear. 

The  new  material  model  with  material  properties  defined  in  Eq.  (5)  was 
incorporated  in  the  SAAS  III  program.  Therefore,  for  an  orthotropic  material, 
the  constants  A^ , B-j,  Ci  and  U0^  are  input  to  the  program  for  the  seven  inde- 
pendent material  properties.  If  the  degree  of  anisotropy  of  the  material  is 
less  than  orthotropic,  e.g.,  transversely  isotropic  or  isotropic,  fewer  con- 
stants (5  sets  for  transverse  isotropy  and  2 sets  for  isotropy)  are  nec  ssary. 

In  SAAS  III,  a maximum  of  six  material  types  are  permitted.  Each  material  can 
have  properties  specified  at  a maximum  of  twelve  temperatures.  Thus,  the 
number  of  sets  of  required  constants  is  a function  of  the  number  of  materials, 
the  anisotropy  of  each  material,  and  the  number  of  temperatures  at  which 
properties  are  specified  for  each  material. 
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2.2  EFFECT  OF  MATERIAL  MODEL  CONSTANTS 

The  general  shapes  of  the  material  property-versus-J  curves  for  different 
.alues  of  the  constants  Ait  , and  are  discussed  in  this  section. 

The  material  model  constant  Aj  is  used  to  define  the  initial  slope  (elastic 
value)  of  the  property.  Thus,  the  effect  of  A-j  on  the  shape  of  the  material 
property-versus-U  curve  can  easily  be  appreciated.  For  Young's  moduli,  values 
of  A are  always  positive  to  guarantee  a positive  strain  energy  density.  For 
Poisson's  ratios,  values  of  A can  be  either  positive  or  negative  because  v can 
have  either  sign  (for  isotropic  materials,  -1<v<1/2;  more  complicated  bounds 
are  given  by  Lempriere  [29]  for  orthotropic  materials).  The  influence  of  the 
material  constants  B-j  and  C-j  and  their  allowable  range  of  magnitude  is  less 
obvious.  Values  of  the  material  property  constants  used  for  Young's  moduli 
for  ATJ-S  graphite  and  7075-T6  aluminum  are  given  ir.  Table  2.  In  general, 
the  values  of  the  material  property  constants  B-j  and  each  range  between 
.05  and  .5  for  both  the  Young's  moduli  and  Poisson's  ratios  of  ATJ-S  graphite. 

A short  study  was  made  to  evaluate  the  effects  of  B-j  and  Ci  on  the  shape  of  the 
material  property-versus-strain  energy  density  curve  and  consequently  on  the 
implied  stress-strain  behavior.  Hereafter,  B^  and  C-j  are  abbreviated  as  B and 
C,  respectively. 


Table  2 Material  Property  Constants  for  the  Young's  Moduli 
of  ATJ-S  Graphite  and  7075-T6  Aluminum 


MATERIAL 

A 

B 

C 

Uo 

ATJ-S  Graphite  (Ez) 
7075-T6  Aluminum  (E) 

1.235  x 106  psi 
10.5  x 106  psi 

. 204557 

2.20069  x 10"4 

.293684 

1.20274 

1 psi 
1 psi 

Before  the  detailed  results  of  the  study  of  the  material  property  con- 
stants B and  C are  discussed,  several  general  comments  are  appropriate.  Only 
positive  values  of  C are  considered  to  be  significant  for  describing  material 
property  behavior.  A negative  or  zero  value  of  C would  not  yield  an  "elastic" 
value  of  material  property  when  the  strain  energy  density  approaches  zero. 
Although  negative  values  of  B have  not  t-en  used  for  ATJ-S  graphite,  they  can 
be  used  in  the  description  of  Young's  roduli  for  strain-hardening  materials, 
i.e.,  concave  upward  stress-strain  curves  as  in  Fio,  10.  Negative  values  of  B 
can  also  be  used  in  the  approximation  of  Poisson's  ratios  for  engineering  metals 
such  as  aluminum  and  steel.  For  these  metals,  Poisson's  ratios  increases  as 
the  strain  energy  density  increases. 

The  schematic  behavior  of  material  property-versus-strair.  energy  density 
curves  for  various  values  of  the  material  property  constants  A,  B,  C,  and  UQ 
is  shown  in  Fig.  11.  As  seen  in  Fig.  11  and  demonstrated  from  the  equations, 
the  family  of  curves  is  coincident  at  a strain  energy  density  value  of  UQ 
(U/UQ  = 1).  The  magnitude  of  B determines  the  ordinate  of  this  intersection. 

In  the  schematic,  the  curves  drawn  are  typical  of  those  for  values  0<B<1  (also 
shown  in  Fig.  11  is  the  straight  line  corresponding  to  the  trivial  case  of  B=0 
and  0<C<  ®).  The  material  property  - strain  energy  relationship  is  linear  for 
Ol.  For  positive  values  of  C less  than  unity,  all  curves  are  concave  downward. 
Specific  examples  of  material  property- vers us -strain  energy  density  curves  for 
values  of  B of  .05  and  .5  are  given  in  Figs.  12  and  13,  respectively.  The 
curves  in  both  figures  are  A = 1 and  UQ  = 1.  At  each  B value,  curves  are  given 
for  four  values  of  C:  .1,  .5,  1,  and  1.25.  Although  no  specific  dimensions 

are  associated  with  the  material  property  constants  in  this  study,  several 
comments  concerning  their  usual  dimensions  are  appropriate. 
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FIGURE  10  STRESS-STRAIN  BEHAVIOR  OF  STRAIN- HARDENING  MATERIALS  WITH  NEGATIVE  B 
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FIGURE  11  SCHEMATIC  BEHAVIOR  OF  MATERIAL  PROPERTY  VERSUS  STRAIN  ENERGY  CURVES 
FOR  VARIOUS  VALUES  OF  A,  B,  C,  AND  U 
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FIGURE  12  MATERIAL  PROPERTY  VERSUS  STRAIN  ENERGY  CURVES  FOR  A-l , B=.05,  U =1 
AND  C=.l,  .5,  1,  AND  1.25  0 
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FIGURE  13  MATERIAL  PROPERTY  VERSUS  STRAIN  ENERGY  CURVES 
FOR  A=1 , B=.5,  U =1,  AND  C=1 , .5,  1,  AND  1.25 
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The  dimensions  of  the  material  property  are  reflected  in  the  dimensions 
of  A.  For  example,  for  Young's  moduli,  A could  have  dimensions  of  pounds  per 
square  inch  (psi)  or  Newtons  per  square  meter  (N/m^).  For  Poisson's  ratios, 

A is  dimensionless.  The  constant  U0  is  included  in  Eq.  5 to  establish  a di- 
mensionless ratio  and,  thus,  a dimensionally  correct  equation  independent  of 
the  alue  of  C.  Both  U0  and  U have  the  dimensions  of  energy  density  (stress). 
Constants  B and  C are  always  dimensionless. 

For  a given  direct  modulus  equation,  the  shape  of  the  corresponding 
stress-strain  curve  can  be  determined.  The  results  for  C = .1,  .5,  1 and 
1.25  are  presented  in  Fig.  14  for  B = .05  and  in  Fig.  15  for  B = .5.  These 
curves  are  in  a normalized  form  (A  = 1 and  U0  = 1).  The  curves  in  Figs.  14 
and  15  are  coincident  at  the  point  corresponding  to  a strain  energy  density 
of  U0  (U/U0  =1).  Also  demonstrated  in  the  figures,  especially  in  Fig.  15, 
is  the  tendency  for  the  stress-strain  curves  to  approach  a straight  line  with 
slope  A(l-B)  as  C approaches  zero  from  the  positive  side.  When  C is  large,  the 
implied  stress  may  not  be  monotonically  increasing  as  a function  of  strain  as 
in  Fig.  16(a).  Instead,  the  stress  level  may  decrease  with  increasing  strain 
after  attaining  a maximum  value  as  in  Fig.  16(b)  or  the  curve  for  C = 1.25  in 
Fig.  15.  The  present  material  model  cannot  be  used  to  represent  the  stress- 
strain  curve  in  Fig.  16(b)  if  the  value  of  U is  the  same  at  any  two  (different) 
stress  levels,  i.e.,  the  material  property-strain  energy  density  curve  must  be 
unique.  The  present  material  model  is  therefore  not  suitable  for  large  strain 
applications  where  the  material  property-U  relationship  is  not  unique.  As  will 
be  discussed  later,  large  value  of  C may  also  lead  to  convergence  difficulties 
in  the  material  model . 

The  direct  modulus  case  is  examined  to  evaluate  the  effects  of  the  rela- 
tive magnitudes  of  material  property  constants  B and  C on  the  general  shape  of 
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the  implied  stress-strain  curve.  If  the  material  property  being  described  with 
Eq,  5 is  a secant  modulus,  the  stress-strain  curves  in  Fiq.  17  are  representa- 
tive of  curves  for  two  values  of  B which  differ  by  an  order  of  magnitude,  e.q., 
curves  for  corresponding  values  of  C in  Fiqs.  14  and  15.  The  only  variable  in 
this  schematic  is  B,  i.e.,  the  material  property  constants  A,  C,  and  U0  are 
fixed.  The  points  shown  on  the  curves  correspond  to  the  implied  stress  and 
strain  values  for  the  condition  U = U0.  As  can  be  seen  in  Fig.  17,  increasing 
B leads  to  a stress -strain  curve  with  more  curvature.  Thus,  the  value  of  B is 
associated  with  the  curvature  of  the  stress-strain  curve.  On  the  other  hand, 
increasing  the  value  of  C in  Figs.  14  and  15  leads  to  an  increased  rate  of 
change  of  curvature  for  a fixed  value  of  B.  Thus,  the  value  of  C is  associated 
with  the  rate  of  change  of  curvature  of  the  stress-strain  curve. 

The  stress-strain  curves  in  Fig.  18  are  representative  of  all  combinations 
of  high  and  low  values  of  B and  C.  The  four  combinations  displayed  in  Fig.  18 
can  be  characterized  as  follows: 

(a)  Low  B,  Low  C - small  initial  curvature  with  slow  curvature  growth 

(b)  High  B,  Low  C - larger  initial  curvature  than  (a)  with  slow  curva- 

ture growth 

(c)  Low  B,  High  C - small  initial  curvature  with  rapid  curvature  growth 

(d)  High  B,  High  C-  larger  initial  curvature  than  (c)  with  rapid  curva- 

ture growth 

In  conclusion,  the  initial  value  of  the  material  property  is  reflected  by 
the  positive  material  property  constant  A.  If  the  material  property  being 
approximated  is  a direct  modulus,  the  value  of  B is  associated  with  the  curva- 
ture of  the  stress-strain  curve.  The  value  of  C,  which  is  positive,  is 
associated  with  the  rate  of  change  of  curvature  of  the  stress-strain  curve. 
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FIGURE  18  STRESS-STRAIN  RESPONSE  FOR  VARIOUS  RELATIVE  MAGNITUDES  OF  B AND  C 


2.3  ITERATION  AND  CONVERGENCE  STUDIES 

The  stresses,  strains,  and  material  properties  are  no  longer  constants  but 
are  highly  interdependent  functions  of  U.  Thus,  an  iterative  procedure  was  de- 
vised to  simultaneously  satisfy,  within  some  small  error,  the  nonlinear  stress- 
strain  relations  and  the  material  property  strain  energy  equations.  This  pro- 
cedure is  schematically  shown  in  Fig.  19.  »he  first  step  is  to  determine  the 
material  property-versus-strain  energy  equations  from  available  uniaxial  data. 

In  order  to  start  the  iterative  procedure,  initial  values  for  the  material  prop- 
erties are  necessary.  The  most  logical  choice  is  the  elastic  or  linear  compo- 
nents of  the  material  properties.  The  next  step  is  to  solve  for  the  stresses 
and  strains  using  these  initial  values  of  material  properties  in  Eq.  (5).  Next, 
the  strain  energy, U,  is  found,  and  new  material  properties  are  determined  for 
that  value  of  U.  Then,  the  stresses  and  strains  are  reevaluated  for  the  revised 
material  properties.  The  material  properties,  stresses  and  strains  are  modified 
until  some  convergence  criterion  for  the  relative  change  in  the  strain  energy  is 
satisfied. 

A short  computer  program  was  developed  to  assess  the  feasibility  of  the 
iterative  approach.  Some  results  obtained  in  the  uniaxial  loading  of  a one- 
element  plane  stress  model  are  shown  in  Fig.  20.  The  element  in  Fig.  21  was 
loaded  with  a stress  of  oz  of  3513  psi.  The  predictions  for  the  first  five 
iterations  of  the  procedure  are  shown  in  Fig.  20.  For  this  problem,  the  proce- 
dure converges  rapidly,  a necessity  for  a procedure  that  is  an  integral  part 
of  a large  finite  element  program.  In  such  a program,  a problem  with  a large 
number  of  elements  can  use  j substantial  amount  of  computer  time  for  each 
iteration.  Convergence  is  defined  to  occur  for  a multi -element  problem  only 
when  all  elements  simultaneously  satisfy  the  convergence  criterion. 

Because  the  present  procedure  does  not  contain  any  criterion  on  which  to 
base  failure  of  a material,  a material  property  equation  can  be  used  far  be- 


FIGURE  20  CONVERGENCE  OF  ITERATION  PROCEDURE 


FIGURE  21  ONE-ELEMENT  PLANE  STRESS  MODEL 


yond  its  meaningful  range.  Two  definitions  of  "meaningful  range"  exist: 

1.  The  model  representation  of  the  material  property  versus 
strain  energy  curve  might  go  to  zero  as  the  energy  in- 
creases such  as  do  the'curves  for  C>1  in  Fig.  11.  That 
is,  we  may  have  a bad  curve  fit  of  the  material  property - 
versus-strain  energy  relation  as  will  be  discussed  in 
Section  2.4.2.  The  computational  consequence  of  the 
material  property-versus-stra in  energy  curve  touching 

or  going  below  zero  is  a zero  or  negative  element  on 
the  diagonal  of  the  compliance  matrix.  This  possibility 
is  countered  by  users  of  the  model  by  the  methods  dis- 
cussed in  Section  2.4.2. 

2.  The  predicted  energy  may  be  greater  than  that  required  to 
fracture  the  material.  We  have  not  included  a failure 
criterion  as  a termination  of  the  present  procedure. 

That  is,  we  predict  strains  wr:ich  could  be  above  failure 
and  hence  beyond  the  meaningful  range.  However,  the 
natural  use  of  the  material  model  in  a design  environ- 
ment is  in  conjunction  with  a failure  criterion.  That 
coupling  of  procedures  will  not  be  discussed  in  this 
report . 

This  effect  of  the  magnitude  of  the  material  constants  on  the  convergence 
of  the  iterative  procedure  was  studied  by  use  of  a brief  computer  program  for 
a "one-element"  uniform  stress  model.  Convergence  of  the  procedure  is  influ- 
enced by  variables  such  as  the  loading,  the  material  properties,  how  the 
material  properties  are  approximated,  and  the  convergence  criterion.  Thus, 
due  to  the  number  of  variables,  the  results  of  this  study  are  somewhat  limit- 
ed. For  this  investigation,  the  material  is  considered  to  be  elastically 
isotropic.  A constant  Poisson’s  ratio  of  .32  is  used  along  with  a constant 
shear  modulus  of  4 x 106  psi . The  convergence  criterion  is  fixed  to  yield 
three  signficant  figure  accuracy  for  the  strain  energy  density.  The  value 
of  the  nonlinear  Young's  modulus  (£)  is  varied  to  evaluate  its  effect  on 
the  convergence  of  the  iterative  procedure.  The  material  property  constants 
A and  U0  are  fixed  at  10.5  x 10&  psi  and  1 psi,  respectively.  The  modulus 
variations  are  induced  by  varying  tne  material  property  constants  E and  C 
as  shown  in  Table  3.  For  each  case,  a uniaxial  tension,  biaxial  tension. 
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and  triaxial  tension  loading  condition  with  equal  normal  stresses  of  72  ksi 
is  investigated. 


For  each  set  of  B and  C values  for  the  three  loading  conditions,  the 
following  results  are  presented  ir  Table  3;  an  indication  of  the  convergence 
of  the  procedure,  the  number  of  iterations  required  for  convergence,  and  the 
strain  energy  density  at  convergence.  For  the  divergent  cases,  the  number 
of  iterations  is  the  number  before  a change  in  sign  of  the  modulus  occurred 
(a  definite  indication  of  divergence).  Obviously,  a strain  energy  density 
value  for  a divergent  case  would  not  be  meaningful. 


Table  3 EFFECT  OF  LOADING  AND  MATERIAL  CONSTANTS  B AND  C ON  CONVERGENCE  OF  ITERATION  PROCEDURE 


CASE 

B 

c 

LOADING 

UNIAXIAL 

BIAXIAL 

TRIAXIAL 

CONVERGE? 

Hi 

U 

CONVERGE? 

NUMBER  OF 
ITERATIONS 

U 

CONVERGE? 

@111 

U 

1 

.001 

31 

.5 

Yes 

2 

247 

Yes 

2 

337 

Yes 

2 

267 

2 

1. 

Yes 

3 

260 

Yes 

4 

362 

Yes 

3 

283 

3 

1.15 

Yes 

4 

290 

Yes 

6 

445 

Yes 

5 

321 

4 

1.2 

Yes 

6 

322 

No 

8 

- 

Yes 

7 

367 

5 

1.25 

No 

10 

- 

No 

3 

- 

No 

6 

- 

6 

.00* 

9001 

1. 

Yes 

2 

247 

Yes 

2 

336 

Yes 

2 

267 

7 

1.8 

Yes 

4 

269 

Yes 

5 

405 

Yes 

4 

295 

8 

1.85 

Yes 

5 

281 

No 

12 

- 

Yes 

5 

313 

9 

1.9 

Yes 

6 

308 

No 

4 

- 

Yes 

10 

372 

10 

1 

f 

1.93 

No 

7 

- 

No 

3 

- 

No 

S 

- 

Three  general  conclusions  can  be  reached  from  the  parametric  study: 

1.  If  only  C changes  in  the  material  property  equation  for  Young's 
modulus,  the  greater  the  value  of  C,  the  slower  the  convergence 
and  the  greater  the  likelihood  of  divergence. 

2.  Similarly,  if  only  B changes  in  the  material  property  eauation 
for  Young's  modulus,  the  greater  the  value  of  B,  the  slower  the 
convergence  and  the  greater  the  likelihood  of  divergence. 

3.  If  only  the  loading  conditions  vary,  the  uniaxial  tension  case  is 

most  likely  to  converge  and  will  converge  most  rapidly  of  the  three 
cases  (uniaxial,  biaxial,  and  triaxial)  examined.  At  the  same  time, 
the  biaxial  tension  case  is  the  least  likely  to  converae  or  will 
converge  the  slowest.  Thus,  biaxial  stress  states  are  the  most 
crucial  test  of  the  convergence  of  the  material  model. 

Thu*;,  from  conclusions  1 and  2,  the  greater  the  nonlinearity  and  the  more 

rapidly  it  grows,  the  greater  the  possibility  of  divergence.  The  third  con- 
clusion may  not  seem  palatable  at  first  glance.  Instead,  the  triaxial  loadinq 
condition  might  be  suspected  to  represent  the  most  severe  test  of  convergence. 
However,  analysis  of  the  energy  conditions  in  uniaxial,  biaxial,  and  triaxial 
loading  states  leads  to  a rational  explanation  of  the  situation. 

Consider  a linear  elastic  body  composed  of  an  isotropic  material  with  a 
Poisson's  ratio  of  .32  and  loaded  uniaxially,  biaxially,  and  triaxially  at  the 
same  stress  level.  The  strain  energy  density  is  .50ce  for  the  uniaxial  stress 
state,  .68v£  for  the  biaxial  stress  state,  and  .54oe  for  tie  triaxial  stress 
state.  Thus,  36%  greater  strain  energy  density  exists  in  the  biaxial  stress 
state  than  in  the  uniaxial  stress  state.  Also,  8%  greater  strain  energy  density 
exists  in  the  triaxial  stress  state  than  in  the  uniaxial  stress  state.  These 

energy  differences  are  larger  for  nonlinear  elastic  bodies  because  the  higher 

the  energy  level,  the  greater  the  strain  nonlinearity.  However,  the  biaxial 
stress  state  obviously  has  the  largest  strain  energy  density.  Therefore,  it 
should  not  be  surprising  that  the  biaxial  stress  state  converges  slower  and 
is  more  likely  to  diverge  than  the  triaxial  or  the  uniaxial  stress  states. 
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2.4  MATERIAL  MODEL  SENSITIVITY  STUDIES 

Three  distinct  sensitivity  studies  are  discussed.  First,  the  symmetry 
or  nonsynmetry  of  the  compliance  matrix  is  investigated.  Second,  different 
approaches  are  examined  for  finding  the  constants  in  the  material  property 
equations  [Eq.  (5)].  Finally,  the  deviations  in  predicted  strains  resulting 
from  +10%  variations  in  uniaxial  strain  data  are  determined. 

2.4.1  Symmetry  of  the  Compliance  Matrix 

In  the  formulation  of  the  material  model , deformations  are  assumed  not  to 
alter  the  degree  or  direction  of  the  anisotropy  of  the  material.  This  assump- 
tion is  reflected  in  the  symmetry  of  the  compliance  matrix  and  its  similar 
form  at  all  deformation  levels.  Since  the  material  model  can  be  generalized 
to  different  degrees  of  anisotropy,  some  of  the  implications  of  this  assump- 
tion for  ATJ-S  graphite  can  be  evaluated. 

The  computer  program  used  to  assess  the  feasibility  of  the  iterative 
approach  was  used  in  this  study.  In  the  program,  a one-element  plane  stress 
model  is  subjected  to  constant  normal  stresses  as  in  Fig.  21.  No  displacement 
boundary  conditions  are  imposed  on  the  model . The  stress-strain  relations  for 
the  case  under  consideration  degenerate  to 


Material  property  versus  U equations  for  the  four  material  properties 
(Ez,  Ee , vz6,  and  vez)  in  Eq.  (6)  were  determined  from  Jortner's  tension  tests 
on  uniaxial  speciments  [30].  The  numerical  values  of  the  constants  in  the 
material  property  equations  used  in  this  study  are  shown  in  Table  4. 


39 


TABLE  4:  BASELINE  VALUES  OF  CONSTANTS  IN  MATERIAL  PROPERTY  EQUATIONS* 
FOR  NONSYWETRIC  COMPLIANCE  MATRIX  STUDY 


PROPERTY 


1. 370x1 06 
1. 720x1 06 
.1607 
.2018 


* A,  B,  C,  and  UQ  are  defined  in  Eq.  (5). 


B 

C 

u0 

.2318 

.2551 

1. 

.1016 

.4458 

1. 

.4456 

.1578 

1. 

.4401 

.1470 

1. 

TABLE  5:  STRAIN  RESULTS  FOR  NONSYMMETRIC  COMPLIANCE  MATRIX  STUDY 


1 1000 
2 2000 

3 3000 

4 4000 

5 1000 

6 2000 

7 3000 

8 4000 

9 0 
10  0 
11  0 
12  0 


Un symmetric 

Symmetric 
with  -vzp/Ez 

Syrnmetric 
wi  "th  —v  a 

cz 

Cc 

Vr 

f 

’-z 

e e 

.000900 

-.000088 

.000900 

-.000088 

.000900 

-.000077 

.002019 

-.000163 

.002019 

-.000163 

.002019 

-.000140 

.003371 

-.000230 

.003371 

-.000230 

.003371 

-.000197 

.005007 

-.000289 

.005007 

-.000289 

.005007 

-.000251 

.000851 

.000550 

.000839 

.000550 

.000852 

.000562 

.001987 

.001244 

.001962 

.001234 

.001989 

.001269 

.003450 

.002143 

.003416 

.002140 

.003455 

.002176 

.005382 

.003385 

.005348 

.003381 

.005388 

.003417 

-.000079 

.000619 

-.000089 

.000619 

-.000079 

.000619 

-.000144 

.001313 

-.000167 

.001313 

-.000144 

.001313 

-.000203 

.002100 

-.000238 

.002100 

-.000203 

.002100 

-.000259 

.003003 

-.000303 

.003003 

-.000259 

.003003 

Three  different  forms  of  the  compliance  matrix. 
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were  used  in  strain  predictions  for  twelve  different  loading  conditions. 

Strain  predictions  for  each  of  the  cases  analyzed  are  given  in  Table  5. 

The  nonsymmetry  of  the  compliance  matrix  about  the  main  diagonal  has  an  almost 
negligible  affect  on  the  strain  predictions  for  the  present  cases.  Thus,  the 
increased  complexity  created  by  the  assumption  of  an  unsymmstric  compliance 
matrix  may  not  be  justified  for  ATO-S  graphite. 
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2.4.2  Evaluation  of  Material  Constants 

The  second  sensitivity  study  was  performed  to  establish  guidelines  for 
evaluation  of  the  cor.otants  appearing  in  the  material  property  equations.  A 
typical  material  property  versus  U curve  for  a direct  modulus  is  shown  as  the 
solid  curve  in  Fig.  22.  The  open  circles  are  the  actual  data  points  used  to 
define  the  curve.  Beside  each  circle  is  the  uniaxial  stress  level  for  that 
point.  Also  shown  are  three  different  numerical  fits  of  Eq.  (5)  to  the  uni- 
axial test  data.  These  three  fits  were  obtained  by  simultaneously  solving  for 
the  constants  B and  C of  Eq.  (5)  using  different  data  points.  For  all  three 
cases,  the  same  initial  value  was  used  for  the  constant  A and  a value  of  unity 
was  assigned  to  !J0.  The  data  points  used  to  define  each  curve  are 
Type  of  Fi t Data  Points  Used 

Low  0,  2000,  3000 

Standard  0,  3000,  5000 

High  0,  4000,  5000 

Three  sets  of  material  property  equations  for  Ez,  E0 , and  vZ0  were  found 

in  the  manner  just  outlined.  The  stress  levels  of  the  data  points  were  not 
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FIGURE  22  APPROXIMATION  OF  M0DULUS-EN3RGY  DATA  WITH  MATERIAL  PROPERTY  EQUATIONS 
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identic,  for  all  properties.  Each  set  of  material  property  equations  was 
analyzed  at  a number  of  loading  conditions  using  the  computer  program  for  the 
one-element  model  previously  described.  The  compliance  matrix  was  assumed  to 
be  symmetric  with  the  terms  off  the  main  diagonal  being  -vze/Ez.  The  constants 
in  the  material  property  equations  were  determined  from  Jortner's  experimental 
results  for  the  :niaxial  test  specimens  [30]. 

Strain  results  are  presented  in  Table  6 in  the  form  of  percentage  devia- 
tions from  strains  obtained  with  the  so-called  standard  fit  of  the  material 
property  equations  to  the  data.  These  percentage  deviations,  which  are  a mea- 
sure of  error,  are  the  difference  between  a ray  from  the  origin  of  a biaxial 
strain  plot  as  in  Fig.  23  to  the  predicted  values  and  a ray  from  the  origin  to 
the  baseline  values.  The  differences  between  the  three  types  of  material 
property  equation  fits  are  more  pronounced  as  the  energy  level  increases.  In 
fact,  for  loading  condition  8,  the  procedure  is  divergent  for  the  low  fit  (this 
procedure  apparently  can  diverge  if  the  material  properties  change  rapidly  with 
U).  The  correlation  between  the  three  r;-Fferent  data  fits  is  relatively  good 
at  low  energy  levels.  Perhaps  the  best  approach  to  establish  the  constants  in 
the  material  property  equations  is  to  use  the  standard  fit  or  high  fit  approach 
in  which  the  approximate  energy  for  convergence  is  determined.  Then,  the  con- 
stants in  material  property  equations  are  adjusted  to  betier  fit  the  material 
properties  in  that  region. 


2.4.3  Variation  in  Material  Properties 

The  third  sensitivity  study  was  undertaken  to  establish  the  deviations  in 
predicted  strain  results  introducted  by  +10%  variations  in  the  input  uniaxial 
stress-strain  data  [30].  At  all  stress  levels,  the  strains  are  incremented  by 
10%  as  shown  in  Fig.  24.  For  the  Poisson's  ratio,  only  the  transverse  strains 
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TABLE  6:  DEVIATIONS  IN  STRAINS  FROM  DIFFERENT  PROCEDURES  FOR  DETERMINING 
MATERIAL  CONSTANTS 


°e 

PERCENT  DEVIATION 
ROM  STANDARD  FIT 

psi 

HIGH  FIT 

LOW  FIT 

0 

1.48 

- 4.41 

0 

1.34 

0.05 

0 

0.81 

10.65 

0 

0.02 

46.35 

1000 

3.09 

- 3.90 

2000 

2.20 

2.82 

3000 

0.57 

27.01 

4000 

-2.35 

* 

1000 

2.09 

- 2.71 

2000 

2.53 

- 2.48 

3000 

2.25 

- 0.03 

4000 

1.49 

6.73 

Procedure  diverges  for  this  case. 


DEVIATION  ■ ^UNFWUE  * 100 
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FIGURE  23  INTERPRETATION  OF  PERCENTAGE  DEVIATIONS  IN  STRAINS 
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are  altered  from  the  baseline  values.  New  constants  are  determined  for  the 
material  property  equations  from  the  artificially  incremented  strains.  Note 
that  all  properties  are  simultaneously  incremented. 

Predicted  strain  results  obtained  using  the  program  for  the  one-element 
model  are  given  in  Table  7.  The  loading  conditions  and  compliance  matrix  sym- 
metry are  that  used  in  the  Evaluation  of  Material  Constants  study.  The  per- 
centage deviations  are  interpreted  as  before  (see  Fig.  23).  A 10%  variation  in 
the  uniaxial  strain  data  used  to  define  the  material  property  eouations  leads 
to  approximately  the  same  numerical  variation  in  the  predicted  strain  results. 
Thus,  the  material  model  is  obviously  very  stable. 


TABLE  7:  DEVIATIONS  IN  STRAINS  FOR  +10%  VARIATIONS  IN  MATERIAL  PROPERTIES 


CASE 

0 . 

z 

e 

psi 

psi 

PERCENT  DEVIATIONS 
FROM  STANDARD  FIT 


0.09 

0.07 

10.04 

10.03 

8.77 

9.12 

9.25 

9.45 

10.07 
10.20 
10.12 

10.08 


-10.07 

-10.05 

-10.03 

-10.01 

- 8.96 

- 9.18 

- 9.35 

- 9.53 
-10.17 
-10.10 
-10.09 
-10.06 
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3.  ELASTIC  MODELS  FOR  DIFFERENT  MODULI 


IN  TENSION  AND  COMPRESSION 


3.1  INTRODUCTION 


The  elastic  behavior  of  materials  with  different  moduli  under  tensile 
loading  than  under  compressive  loading  is  described  briefly  in  Section  1.2. 
Two  material  models  for  this  behavior  will  be  described  in  this  section. 

In  both  models,  the  compliances,  S^j,  in  the  strain-stress  relations  for 
axi symmetric  bodies  under  axi symmetric  loading 
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are  determined  based  on  the  principal  stress  state.  In  principal  stress  co- 
ordinates, the  strain  stress  relations  are 
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Note  that  (1)  the  e-direction  is  always  a principal  stress  direction  because 
of  the  axisymmetric  loading  and  (2)  principal  stress  directions  and  principal 
strain  coordinates  do  not  coincide  for  orthotropi'  ..terials  as  they  do  for 
isotropic  materials.  The  relation  between  the  principal  stress  (p-q-e)  coor- 
dinates and  the  body  (r-z-9)  coordinates  is  shown  in  Fig.  25.  The  e-direction 
is  toward  the  reader  out  of  Fig.  25.  Note  that,  by  definition,  the  body  coor- 
dinates coincide  with  the  principal  material  coordinates  for  graphite.  Trans- 
formation of  stresses,  strains,  and  material  properties  between  these  coordi- 
nate systems  will  be  necessary. 

The  objective  in  both  material  models  is  to  define  a rational  procedure 
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FIGURE  25  RELATION  OF  PRINCIPAL  STRESS  (p-q 
TO  BODY  AND  PRINCIPAL  MATERIAL  (r 
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for  assigning  the  S^j  on  the  basis  of  the  principal  stress  state  and  the  values 


of  the  S??  under  tensile  loading  and  those  values  under  compressive  loading. 


Because  the  material  properties  depend  on  the  stress  state  and  vice  versa,  the 
basic  problem  is  statically  indeterminate.  However,  the  indetermi nancy  can  be 
resolved  by  an  apparently  convergent  iterative  procedure  consisting  of  four 
steps.  First,  displacement  and  stress  calculations  are  performed  based  on  an 
initial  assumption  of  stress  signs  with  an  implied  initial  choice  of  material 
properties.  Second,  the  appropriate  material  properties  are  selected  based  on 
the  principal  stresses  calculated  in  the  previous  step.  Third,  displacements 
and  stresses  including  the  new  principal  stresses  are  recalculated.  Fourth, 
steps  two  and  three  are  repeated  until  convergence  to  tne  desired  accuracy  is 
achieved. 

The  first  model,  called  the  Weighted  Compliance  Model,  is  due  to  Jones 
[22,  24].  This  model  consists  basically  of  adding  the  tensile  ’.nd  compressive 
compliances  in  proportion  to  the  presence  of  the  respective  tensile  and  com- 
pressive principal  stresses;  hence,  the  name  Weighted  Compliance  Model  is  used, 
fio  theoretical  basis  exists  for  the  Weighted  Compliance  Model;  it  is  an  engi- 
neering approach,  rather  than  a scientific  approach,  to  a very  difficult  prob- 
lem. 

The  second  model,  called  the  Restricted  Compliance  Model,  is  due  mainly 
to  Isabekian  and  Khachatrian  [23]  (modifications  necessary  to  extend  Isabekian 
and  Khachatrian's  model  from  a plane  stress  state  to  an  axi symmetric  stress 
state  are  described  in  Section  3.3).  This  model  has  a more  scientific  basis 
than  the  Weighted  Compliance  Model.  Specifically,  the  compliance  matrix  is 
made  symmetric  by  prescribing  certain  relations  between  the  tensile  and  com- 
pressive properties  so  that  they  satisfy  the  known  transformations  of  aniso- 
tropic elasticity.  However,  by  enforcing  these  relations  between  tensile  and 
compressive  properties,  we  limit  our  ability  to  treat  real  engineering  materials, 
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That  is,  some  materials  with  different  moduli  in  tension  and  compression  do  not 
nave  properties  that  satisfy  the  Isabekian  and  Khachatrian  restrictions. 

Both  models  are  attempts  to  treat  materials  with  different  moduli  in  ten- 
sion and  compression  as  special  anisotropic  elastic  materials  with  the  com- 
pliance matrix  symmetry  characteristics  of  ordinary  elastic  materials  (for 
which  a potential  function  exists).  This  compliance  matrix  symmetry  is  a 
correction  to  the  basic  Ambartsumyan  model  with  a nonsynmetric  compliance 
matrix  [21].  This  correction  is  necessary  if  anisotropic  elasticity  theory 
is  to  be  used  once  the  compliances  are  determined.  Both  material  models  are 
described  in  the  following  sections  and  are  contrasted  in  Section  3.4. 


3.2  WEIGHTED  COMPLIANCE  MODEL 


For  an  axi symmetric  body  under  axi symmetric  load,  the  following  strain- 
stress  relations  in  principal  stress  (p-q-e)  coordinates  were  proposed  by  Jones 
[24]  for  orthotropic  materials  that  exhibit  different  moduli  in  tension  and 


compression: 


rpq  c-pq 


iLsra  ^ si3 

l spq  spq  spq  spq 

•-  13  23  33  36 


Note  again  that  the  principal  stress  directions  do  not  coincide  with  the 
principal  strain  directions.  The  compliances,  S^j,  are  assigned  according  to 

the  signs  and  magnitudes  of  the  principal  stresses: 

if  - , n.  c-pq . cpq 
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if  ''a  > s§!S  = s§§t 
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The  weighting  factors,  kppq,  etc.,  are  used  to  make  the  compliance  matrix 
symmetric.  Some  other  function  of  the  principal  stresses  could  be  used  to  de- 
fine the  weighting  factors.  Full  qualification  of  the  form  of  the  weighting 

factors  awaits  definitive  experimental  work.  Note  that  only  two  of  the  three 

pq  pq 

principal  stresses  are  used  to  determine  each  of  the  cross-compliances  S^.  S^, 
pq 

and  S23.  Furthermore,  a single  principal  stress  is  used  to  determine  each  of 
the  cross-compliances  S^g,  sjjjj,  and  S^g.  The  compliance  Sgg  cannot  be  ration- 
ally defined  nor  is  necessary  for  transformation  to  any  other  coordinate 

pq 

system.  The  reason  for  the  perhaps  surprising  lack  of  importance  of  S66  1S 

that  it  vanishes  identically  from  all  transformation  relations  from  principal 

stress  coordinates  to  any  other  coordinates. 

The  compliance  and  sP*?  in  the  principal  stress  (p-o)  coordinates 
1 J c * jc 

arc  related  to  tne  compliances  Srz  and  Srz  in  the  principal  material  and 

Ut  ijc 

body  (r-z)  coordinates  by  the  usual  transformations  of  anisotropic  elasticity: 


2.  , erz 
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■ sntc>  5i"6COSB 


(12) 

(cont'd.) 


where  the  subscript  t or  c is  taken  as  appropriate,  and  8 is  the  angle  be- 
tween the  body  (r-z-e)  and  principal  material  coordinates  and  the  principal 
stress  (p-q-e)  coordinates  as  defined  in  Fig.  25. 

rz 

The  compliances  in  principal  material  coordinates,  S-.^c,  are  related 
to  the  engineering  constants  (direct  moduli,  Poisson's  ratios,  and  shear 
moduli)  by: 


Vztc^rtc  - -vzrtc^ztc 


v*  = - /F 
13tc  Vetc  rtc 

S22tc  ’ 1/Eztc 


’vertc^etc 


Sr^  = /£ 

23tc  'zetc  ztc 


vsztc^6tc 


S33tc  = 1/Eetc 


S',  = 1/G  . 

66tc  rztc 


where  . . - -l  /<  for  = c.  and  all  other  stresses  are  zero.  There  are 

rzt  z r r t 

apparently  seven  independent  material  properties  in  tension  in  Eq.  (13)  and 

rz  rz 

the  same  number  in  compression.  However,  the  compliances  S and  Sggc 
d/Grzt  and  VGrzc, respectively)  cannot  be  measured  in  a shear  test  on  an 
orthotropic  material  with  different  moduli  in  tension  and  compression  since 
one  principal  stress  is  tension  and  the  other  is  compression.  Instead,  in 
accordance  with  a suggestion  by  Tsai  [31],  the  tension  modulus  at  45  degrees 
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to  the  principal  material  axes,  E^,  is  measured  and  Sg|t  is  obtained  from 

2v 


5rz  = _L 
66t  T 
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El 


r t 


1-  + _2 
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rzt 


(14) 


zt 
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rt 


A similar  relation  is  used  to  define  Sg*c  in  terms  of  E^. 

For  axi symmetrically  loaded  axi synwetri c transversely  isotropic  bodies, 


Ertc  = £6tc  vrztc  = veztc  S0  remain^n9  independent  properties  are 


Ertc*  vrztc*  Vetc’  ^ztc*  ^tc 


(15) 


i.e.,  five  properties  in  tension  and  five  in  compression  for  a total  of  ten 
independent  properties. 

3.3  RESTRICTED  COMPLIANCE  MODEL 

Isabekian  and  Khachatrian  [23]  define  a synmetric  compliance  matrix  for 
orthotropic  materials  with  different  moduli  in  tension  and  conpression.  They 
require  the  material  properties  in  principal  material  directions  to  have  inter- 
relationships such  that  the  compliances  are  symmetric  in  any  coordinate  system. 
These  relationships  are  obtained  in  addition  to  the  usual  Ambartsumyan  approach 
of  superposing  stress  states  of  uniaxial  tension  in  one  direction  and  uniaxial 
compression  in  the  other  direction. 

For  example,  if  cp>0  and  cq<0,  then  the  strain-stress  relations  are 

f = spq  , + spq  . *\ 

p 'p  12c  "q 


= spq  _ + cpq  0 
•q  12t  'P  22c  °q 


= sP^  - + S'  ■ c 

!pq  let  Jp  26c°q  ^ 


:pq 


(16) 


For  these  relations  to  be  symmetric. 


spq  = spq 

12t  12c 


(17) 
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if  trie  compliance  matrix  is  to  be  symmetric  in  principal  material  (r-z)  direc- 
tions, then 


"rzt  _ rzc  ’zrt  _ v'zrc 

' £ ~ E 

rc  zt  2c 


Tt  "rr  “zt 

If,  moreover,  £q.  (17)  is  to  Le  valid  in  all  coordinate  system,  then 


(18) 


1 _1_ 
E„*  ' E, 
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i _ 1 


09) 


-rt  fcrc  £zt  Lzc  E|  Ec 
in  whicn  E^  and  Ec  are  trie  tensile  and  compressive  moduli,  respectively,  at  45° 

i I 

to  the  principal  material  directions.  Note  that  Et  and  Ec  are  related  to  the 
shear  moduli  in  all  tension  or  all  compression  stress  states  by  Eq.  (14)  and  its 
compression  equivalent.  Accordingly,  the  five  independent  material  properties 
are 

Ert‘  Ezt’  'rzt’  F't»  Erc  ^ 

or 

Erc*  Ezc’  'rzc'  Ec!  Ert  ^21* 


or  other  appropriate  combinations  of  the  foregoing  relations. 

For  axisymmetric  bodies  under  axisymmetric  loads,  the  foregoing  properties 
for  the  plane  stress  case  must  be  supplemented  by  additional  material  properties. 
Because  of  the  two  kinds  of  axisynmetry,  no  additional  shear  modulus  - related 
properties  are  necessary.  However,  the  direct  modul  \ Eetc  must  be  defined  in 
addition  to  the  Poisson's  ratios  *r9tc  and  v2a^c-  ’The  s-direction  is  automati- 
cally a principal  stress  direction.  Thus,  033  in  any  coordinate  system  is  af- 
fected only  by  E . and  no  other  moduli  because  of  the  simplicity  of  the  Irans- 
J v stc 

formation  relation.  The  cross-compliances  S13  and  S23  are  symmetric  under  mixed 
tensile  and  compressive  loading  in  principal  material  coordinates  if 


WEzt  = Vrec/Ezc 


(22) 
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^zot  = vzec  (23) 

Ezt  Ezc 

However,  Ert>  £rc,  E2t>  and  £2C  are  either  independent  properties  or  dependent 
and  determined  frori  Eq.  (19).  Thus,  only  two  Poisson's  ratios  are  irdependent, 
say  . r t and  . t.  All  transformations  of  stresses  to  obtain  principal  stresses 
are  rotations  about  the  G-direction.  Hence,  no  conditions  are  imposed  cn  E0t 
and  £>c.  That  is,  no  relation  like  Eq.  (19)  exists  between  F0t  and  E0C.  Accord- 
nine  independent  material  properties  are  required  for  analyses  of  axi- 
syimnet.  rally  loaded  axisymmetric  solids  made  of  orthotropic  multimodulus  mate- 
rials. Five  of  these  properties  come  from  the  plane  stress  model  and  four  are 
from  tne  extension  to  axisymmetric  solids.  For  transversely  isotropic  materials 
like  ATJ-S  graphite  with  the  plane  of  isotropy  in  the  r-e  plane,  Ee=Er  and 
rz=-yz»  etCi  With  these  reductions,  the  six  independent  material  properties 
are 


Ert’  Ezt’ 


'rzt’ 


’rot’  Et>  Erc 


(24) 


or  appropriate  equivalents  through  the  foregoing  relations. 

As  one  consequence  of  the  Restricted  Compliance  Model,  the  shear  modulus 
for  pure  shear  in  coordinates  45v  from  the  principal  material  coordinates  can  be 
shown  to  be  independent  of  the  gn  of  the  shear  stress.  Consider,  for  example, 
the  elements  from  a uni  directionally  reinforced  composite  material  shown  in  Fig. 
26.  There,  the  positive  shear  stress  ^ads  to  tensile  principal  stress  in  the 
fiber  direction  and  compressive  prir  jal  stress  in  the  direction  transverse  to 
the  fibers.  Similarly,  the  neg:.f.'  hear  stress  leads  to  compressive  stress  in 
tne  fiber  direction  and  tensile  stress  in  the  transverse  direction.  Because  the 
tensile  ano  compressive  moduli  are  different  in  the  fiber  and  transverse  direc- 
tions, the  shear  moduli  at  45^  would  seem  to  be  different  on  purely  rational 
grou-.ds.  however,  because  of  the  relations  between  tensile  and  compressive 
moduli  in  Eq.  (19),  the  expressions  for  the  shear  modulus  at  45c  under  positive 
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POSITIVE  SHEAR  STRESS 


NEGATIVE  SHEAR  STRESS 


FIGURE  26  POSITIVE  AND  NEGATIVE  SHEAR  STRESS 

AT  45°  TO  PRINCIPAL  MATERIAL  DIRECTIONS 
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shear  stress  can  be  shown  to  be  identical  to  the  expression  for  the  shear 
modulus  at  45°  under  negative  shear  stress. 


3.4  SUMMARY 

Irrespective  of  the  material  model  used,  the  strain-stress  relations  in 
principal  stress  coordinates  must  be  transformed  to  the  body  (r-z-t)  coordinates 
for  solution  of  equilibrium  problems.  Those  transformations  take  place  accord- 
ing to  the  usual  transfon.-.ati on  relations  of  anisotropic  elasticity: 
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Moreover,  stresses  and  displacements  from  both  mooels  for  two  distinct 
problems  cannot  be  superimposed  as  can  solutions  for  linear  elastic  problems. 
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Superposition  is  invalid  because  the  principal  stress  directions  for  two  sets 
of  stresses  on  the  same  body  are  generally  different.  However,  if  the  principal 
stress  directions  are  identical  in  both  solutions,  then  the  solutions  can  be 
superimposed.  Ambartsumyan  [21]  terms  this  general  invalidity  of  superposition 
a form  of  nonlineari ty.  In  practical  terms,  each  problem  must  be  separately 
•rived,  and  no  building  up  of  many  solutions  by  superposition  of  simple  solu- 
tions should  be  expected. 

The  two  mater;-,  models  described  in  this  section  have  different  numbers 
of  independent  material  properties.  The  Restricted  Compliance  Model  has  six 
independent  properties  for  an  axi symmetrically  loaded  axi symmetric  transversely 
isotropic  body.  For  the  same  body,  the  Weighted  Compliance  Model  has  ten  in- 
dependent properties.  Both  models  have  advantages  (and  disadvantages).  The 
Restricted  Compliance  Model  is  derivable  on  a rational  scientific  basis;  however, 
the  compliance  restrictions,  as  will  be  seen  in  Section  4,  are  too  stringent 
for  many  real  materials.  On  the  other  hand,  the  Weighted  Compliance  Model  is 
not  derivable  because  the  weighting  factors  are  somewhat  arbitrary;  however, 
this  model  has  the  flexibility  in  "acceptable"  compliances  to  be  applicable  to 
many  real  materials.  Thus,  the  Weighted  Compliance  Model  is  an  engineering 
approximation  whereas  the  Restricted  Compliance  Model  is  a more  scientific  de- 
scription of  the  basic  phenomenon  , but  is  restricted  to  a more  limited  class 
of  materials. 
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4.  NONLINEAR  DEFORMATION  MODELS  FOR 
DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION 

4.1  INTRODUCTION 

The  objective  of  this  section  is  to  simultaneously  model  the  tensile  and 
compressive  stress-strain  behavior  of  graphite  shown  in  Figs.  27  and  28,  re- 
spectively. The  tensile  behavior  is  typical  of  a material  that  softens  under 
increasing  stress.  That  is,  among  other  characteristics,  the  Poisson's  ratios 
decrease  as  the  stress  increases  leading  to  the  softening  behavior  illustrated 
in  Fig.  3.  On  the  other  hand,  the  compressive  behavior  is  more  nearly  charac- 
teristic of  a hardening  material.  That  is,  the  Poisson's  ratios  are  constant 
or  increase  as  the  stress  increases.  Thus,  we  might  expect  the  hardening  be- 
havior (although  certainly  not  hardening  in  the  sense  of  a concave  upward 
stress-strain  curve  as  in  Fig.  10)  illustrated  in  Fig.  3.  Accordingly,  strain 
profiles  for  various  increasing  stress  levels  might  exist  as  schematically 
shown  in  Fig.  29  with  biaxial  tension  in  the  upper  right-hand  quadrant  and  bi- 
axial compression  in  the  lower  left-hand  quadrant.  However,  no  biaxial  com- 
pression data  exist  to  verify  this  speculation.  Moreover,  the  behavior  is 
not  well-defined,  although  some  measurements  have  been  made  and  will  be  dis- 
cussed in  Section  6. 

The  contrasting  tensile  and  compressive  behavior  illustrated  in  Fig.  29 
may  be  perfectly  realistic  when  the  deformation  mechanisms  in  tension  are  com- 
pared with  the  me cn an isms  in  compression.  Specifically,  microcrackirg  in  ten- 
sion leads  to  an  apparent  increase  in  volume  (decreasing  Poisson's  ratios). 
However,  in  compress  ion,  microcracking  may  not  occur.  The  m crocracki no  phenom- 
enon may  be  related  to  the  porous  nature  of  graphite.  Tne  pores  have  a tendency 
to  open  and  perhaps  tear  under  tensile  stress.  On  the  other  hand,  tne  pores 
would  have  a tendency  to  close  and  perhaps  collapse  under  compression.  Because 
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FIGURE  28  COMPRESSIVE  STRESS-STRAIN  BEHAVIOR  OF  GRAPHITE 
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of  tne  pores,  graphite  and  other  porous  nonmetallic  materials  can  have  plastic 
volume  changes  when  subjected  to  tensile  or  compressive  stresses.  Under  com- 
pression, porous  materials  can  compact  (decrease  in  volume)  whereas,  under 
tension^  they  can  dilate  (increase  in  volume).  Thus,  classical  plasticity 
theories  with  the  usual  zero  plastic  volume  cnange  hypothesis  are  not  appli- 
cable to  porous  materials. 

The  nonlinear  material  model  described  in  Section  2 is  combined  in  this 
section  with  the  linear  multimodulus  material  models  of  Section  3 to  obtain  non- 
linear multimodulus  material  models.  With  the  material  models  of  this  section, 
both  the  nonlinear  and  multimodulus  deformation  behavior  of  an  orthotropic  mate- 
rial can  be  modeled  simultaneously,  a necessity  for  description  of  the  behavior 
of  ATJ-S  oraphite. 

Two  nonlinear  multi modulus  models  with  the  same  basic  approach  to  nonlinear 
behavior  are  considered.  Thus,  the  distinction  between  the  models  is  based 
solely  on  the  multimodulus  formulations  or,  more  specifically,  on  the  deter- 
mination of  the  compliance  matrix.  Isabekian  and  Khachatrian's  extension  [23] 
of  Ambartsumyan's  linear  multimodulus  approach  [21]  is  used  in  what  is  called 
the  "nonlinear  restricted  compliance  model"  (RCM).  The  linear  multimodulus 
compliance  approach  of  the  RCM  model  is  discussed  in  Section  3.2.  In  the  other 
material  model,  called  the  "nonlinear  weighted  compliance  model"  (WCM),  the 
linear  multimodulus  approach  developed  by  Jones  [22,  24]  (which  has  been  de- 
scribed in  Section  3.3)  is  combined  with  the  nonlinear  approach  described  in 
Section  2. 

The  discussion  of  the  nonlinear  multi  modulus  material  models  is  divided 
into  five  major  sections.  The  general  characteristics  of  the  iteration  proce- 
dure common  to  both  material  models  are  discussed  in  Section  4.2.  The  energy 
functions  used  to  relate  the  multiaxial  stress  state  to  the  mater. al  properties 
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are  described  in  Section  4.3.  The  distinctive  characteristics  of  the  nonlinear 
restricted  compliance  materia]  model  and  the  nonlinear  weighted  compliance  model 
are  discussed  in  Section  4.4  and  4.5,  respectively.  Some  summary  remarks  on  the 
two  nonlinear  multimodulus  models  are  given  in  Section  4.6. 

4.2  ITERATION  PROCEDURE  FOR  MATERIAL  MODELS 

The  overall  iteration  procedure  common  to  both  the  nonlinear  multimodulus 
material  models  will  be  described.  In  both  the  nonlinear  model  and  the  multi- 
modulus models  previously  discussed,  an  indeterminate  system  is  solved  with  an 
iterative  approacn.  The  stresses  and  strains  are  dependent  on  the  material  pro- 
perties, and  they,  in  turn,  are  dependent  on  the  stresses  and  strains.  In  the 
nonlinear  approach  of  Section  2,  the  material  properties  are  related  to  the 
strain  energy  (density),  which  is  found  from  the  product  of  the  stresses  and 
the  respective  strains.  In  the  linear  multimodulus  models,  the  composition  of 
the  compliance  matrix,  hence  the  material  properties,  is  determined  from  the 
signs  of  the  principal  stresses.  In  the  nonlinear  multimodulus  material  models, 
the  selection  of  material  properties  and,  consequently,  the  stress-strain  re- 
lationships is  based  on  both  the  signs  of  the  principal  stresses  and  the  mag- 
nitude of  an  energy  function. 

The  iteration  procedure  devised  to  simultaneously  satisfy  the  constraints 
of  *-oth  the  nonlinear  and  multimodulus  problems  is  illustrated  schematically  in 
Fig.  30.  Each  step  in  the  procedure  is  described  in  the  following  paragraphs. 

The  first  step  in  the  iteration  procedure  is  to  determine  the  material 
property-versus-strain  energy  relationships  independently  in  tension  and  in  com- 
pression from  available  uniaxial  data.  For  ATJ-S  graphite,  the  general  form 
of  Eq.  (5)  is  suitable  for  representing  the  material  property  variations  with 
strain  energy.  Thus,  sets  of  constants  A,  B,  C and  U0  are  determined  for  each 
independent  material  property.  Separate  sets  of  constants  are  determined 
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for  the  tension  and  compression  representation  of  a material  property  variation 
with  strain  energy  if  both  sets  are  input  to  the  procedure  as  independent  mate- 
rial properties.  The  nunber  of  independent  properties  required  for  an  analysis 
is  dependent  on  the  symmetry  of  the  body  and  loading,  the  material  model  used, 
and,  of  course,  the  degree  of  anisotropy  of  the  material.  For  an  axisymnetric 
body  under  axisymmetric  load,  the  number  of  independent  material  properties  for 
three  levels  of  anisotropy  (i.e.,  isotropic,  transversely  isotropic,  and  ortho- 
tropic) is  defined  in  the  discussions  of  the  multimodulus  models. 

To  initiate  the  iteration  procedure,  the  linear  or  "elastic"  components 
of  the  tensile  material  properties  are  used.  The  desired  stresses  and  strains 
(including  the  principal  stresses  and  strains)  can  be  computed.  The  principal 
stresses  and  their  orientation  are  required  for  the  multimodulus  formulations. 
The  strain  energies  are  calculated  from  the  stresses  and  strains  by  use  of  one 
of  three  strain  energy  approaches  discussed  in  the  next  section.  The  strains 
in  principal  stress  directions  need  be  computed  only  if  the  separate  contribu- 
tions to  the  total  energy  of  the  tensile  and  compressive  stresses  are  used. 

After  the  strain  energies  are  evaluated,  new  independent  material  proper- 
ties are  found.  This  revised  set  of  material  properties  is  used  to  formulate 
new  all-tension  and  all -comp res  si  on  compliance  matrices.  How  these  matrices 
are  found  is  determined  by  which  multimodulus  approach  is  used.  For  the  same 
material,  fewer  material  properties  are  required  in  the  RCM  than  in  the  WCM 
because  relationships  exist  in  the  RCM  between  the  all-tension  and  the  all- 
compression  compliances.  The  compliance  matrices  for  the  tensile  and  compres- 
sive stress  states  are  then  transformed  to  principal  stress  directions.  Next, 
the  composition  of  the  multimodulus  compliance  matrix  is  determined  by  use  of 
the  signs  of  the  principal  stresses.  The  multimodulus  compliance  matrix  is 
obtained  from  the  transformed  tensile  and  compressive  compliance  matrices  by 
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either  the  weighted  compliance  or  restricted  compliance  approach-  Upon  defi- 
nition of  a multimodulus  compliance  matrix,  the  stresses,  strains,  and  strain 
energies  are  reevaluated.  The  loop  from  the  calculation  of  new  material  prop- 
erties to  the  computation  of  stresses,  strains,  and  strain  energies  is  repeated 
until  the  convergence  criterion  is  satisfied. 

Convergence  is  achieved  in  the  linear  multimodulus  model  developed  by 
Jones  [22,  24]  when,  on  successive  iterations,  the  absolute  value  of  the  rela- 
tive change  in  magnitude  of  each  significant  stress  is  less  than  or  equal  to  a 
specified  value.  The  word  "significant"  denotes  that  the  absolute  value  of 
the  stress  is  large  compared  tc  the  summation  of  the  absolute  values  of  all 
the  stresses.  For  example,  consider  the  axisymmetric  loading  of  an  axisym- 
metric  body.  The  stress  cr  is  "significant"  if 
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Although  such  a convergence  criterion  has  been  demonstrated  to  be  suitable  for 
the  linear  multimodulus  procedure,  it  would  be  ineffective  as  a convergence  in- 
dicator of  the  present  nonlinear  procedure.  Consider  the  example  of  Fig.  20 
in  which  the  stress  is  constant  and  only  the  strain  varies  on  successive  iter- 
ations. Convergence  would  be  incorrectly  indicated  on  the  second  iteration 
with  the  stress-based  criterion  since  the  stress  is  unchanged. 

A convergence  criterion  based  on  the  relative  change  in  total  strain 
energy  is  compatible  with  both  a nonlinear  and  a multimodulus  material  model 
and,  therefore,  is  used  in  nonlinear  multimodulus  material  models.  Whenever 
the  absolute  value  of  the  relative  change  in  total  strain  energy,  |REL  aU|, 
between  two  iterations  i-1  and  i 
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(28) 


|RR.  AU|  = l(Ui-U._1)/Ui_1 i 
becomes  sufficiently  small,  then  the  iterative  procedure  terminates  and  conver- 
gence is  assumed. 

When  the  nonlinear  multimodulus  material  models  are  used  in  finite  ele- 
ment computer  programs,  convergence  is  defined  in  this  report  to  occur  when 
the  convergence  criterion  is  simultaneously  satisfied  by  each  element.  A fi- 
nite element  problem  with  a large  number  of  elements  can  require  a substantial 
amount  of  computer  time  for  each  iteration.  Thus,  the  nonlinear  multimodulus 
procedures  were  specifically  devised  to  obtain  rapid  convergence  by  considering 
both  the  nonlinear  and  multimodulus  characters  of  the  problem  on  the  same  iter- 
ation. 

4.3  STRAIN  ENERGY  FUNCTIONS 

For  the  nonlinear  multimoaulus  material  models,  three  different  energy 
functions  are  investigated:  t:tal  strain  energy,  divided  strain  energy,  and 

weighted  strain  energy.  Each  of  these  functions  will  be  discussed  brieily. 

The  strain  energies  used  in  che  material  models  are  actually  strain  energy 
densities.  Although  the  word  density  is  frequently  omitted,  the  strain  ener- 
gies are  always  understood  to  be  on  a unit  volume  basis. 

The  total  strain  energy  defined  in  Eq.  (4)  can  be  used  to  determine  all 
material  properties  in  each  iteration  as  is  done  in  the  nonlinear  model  of 
Section  2.  Thus,  application  of  the  total  energy  function  is  the  same  for  both 
the  nonlinear  material  model  and  the  nonlinear  multimodulus  material  models. 

In  the  divided  energy  approach,  the  total  strain  energy  is  separated  into 
two  components:  (1)  the  contribution  from  the  tensile  principal  stresses  and 

(2)  the  contribution  from  the  compressive  principal  stresses.  These  two  com- 
ponents are  not  invariant  under  rotation  of  coordinates  and  are  defined  only 
in  principal  stress  coordinates.  However,  their  sun,  the  total  strain  energy. 
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is  invariant  witn  respect  to  coordinate  transformations.  In  the  divided  ener- 
gy approach,  the  tensile  component  of  strain  energy  is  used  to  determine  the 
tensile  material  properties,  and  the  compressive  component  of  strain  energy  is 
used  to  determine  the  compressive  material  properties. 

Different  nonlinear  stress-strain  response  in  uniaxial  tension  and  in  uni- 
axial compression  for  ATJ-S  graphite  was  found  by  Jortner  [30]-  The  different 
deformation  behavior  under  tensile  and  compressive  loads  might  be  related,  at 
least  in  part,  to  a microcracking  phenomenon  (to  which  plastic  volitne  changes 
have  also  been  attributed).  In  the  nonlinear  multiraodulus  material  models, 
multiaxial  stress  states  are  plated  to  the  material  properties  defined  in  uni- 
axial stress  states  by  a strain  energy  function.  The  relationship  between 
strain  energy  and  material  properties  might  be  expected  to  be  dependent  on  the 
signs  of  the  multiaxial  principal  stresses.  The  division  of  energy  is  a mech- 
anism used  to  investigate  such  a dependency. 

In  the  weighted  energy  approach,  the  effective  energy  level  in  terms  of 
the  tensile  and  compressive  components  of  the  total  strain  energy  density  is 

uw=(;rj'uc  t[,'(r)”]Ut  (29) 

where 

n = a positive  integer  constant 
U = total  strain  energy 
Uc  = strain  energy  of  compressive  stresses 
Ut  = strain  energy  of  tensiie  stresses 
U.w  = weighted  energy 

This  energy  level,  U^,  is  used  to  find  both  the  ti—iile  and  the  compressive 
material  properties.  Thus,  this  approach  is  more  desirable  than  the  divided 
energy  method  from  the  standpoint  that  a consistent  energy  level  is  used  to 
determine  the  current  values  of  all  material  properties.  When  n=l  in  Eq.  (29), 
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the  weighted  energy  density  varies  between  the  average  and  the  sum  of  Ut  and  Uc 
as  a function  of  the  ratio  Uc/U.  This  variation  is  shown  in  Fig.  31  along  with 
the  variation  for  n=2.  When  n=l,  the  variation  is  symmetrical  about  Uc/U=.5. 

On  the  other  hand,  for  n=2,  the  weighted  energy  density  is  larger  when  ( Uc : ) = 
(1:9)  than  when  (Uc:Ut)  = (9:1).  When  the  weighted  energy  is  used  in  the  com- 
parisons between  predicted  and  measured  strains,  n is  always  unity. 

jU  NONLINEAR  WEIGH : £D  COMPLIANCE  MATERIAL  MODEL 

Three  aspects  of  the  nonlinear  ..eighted  compliance  material  model  (WCM) 
are  discusseu:  (1)  the  distinctive  features  of  the  iteration  approach,  (2) 

the  restrictions  imposed  on  the  stress-strain  behavior  that  can  be  treated 
with  the  model,  and  (3)  the  computerized  application  of  tse  material  model. 

Two  steps  in  the  general  nonlinear  multi  modulus  procedure  outlined  in  Fig. 
30  are  dependent  on  which  multimodulus  procedure  is  used.  The  first  step  is 
the  evaluation  of  the  all-tension  and  all -compress ion  compliance  matrices.  The 
elements  of  these  matrices  are  determined  from  the  material  properties  by  the 
same  expressions  in  both  the  linear  and  nonlinear  multi  modulus  approaches.  In 
the  linear  approach  of  Section  3.2,  the  material  properties  are  constant.  How- 
ever, in  the  nonlinear  approach,  the  material  properties  are  functions  of  the 
strain  energy.  The  second  step  which  depends  on  the  multimodulus  procedure  is 
the  computation  or  the  multimodulus  compliance  matrix.  This  step  is  identical 
in  both  the  linear  and  the  nonlinear  weighted  compliance  material  models. 

When  the  linear  multimcdulus  approach  cf  Section  3.2  is  combined  with  the 
nonlinear  approach  of  Section  2,  restrictions  are  placed  on  the  deformation 
behavior  of  the  material  that  can  be  treated.  As  with  the  nonlinear  material 
model  of  Section  2,  the  reciprocal  relations,  Eq.  (3),  between  material  proper- 
ties '-.re  assumed  to  be  valid  when  all  mater;..l  properties  are  evaluated  at  the 
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same  energy  level.  Reciprocal  relations  exist  for  both  tensile  and  compres- 
sive properties  in  the  multimodulus  problem.  Symmetry  of  the  multimodulus  com- 
pliance matrix  in  each  iteration  is  ensured  in  the  nonlinear  WCM  by  weighting 
factors  used  in  the  development  of  the  multimodulus  compliance  matrix,  not 
by  relationships  between  tensile  and  compressive  material  properties  as  will  be 
done  in  the  RCM.  Thus,  with  the  nonlinear  weighted  compliance  material  model, 
fewer  restrictions  are  placed  on  the  stress-strain  behavior  that  car.  be  treated. 
However,  in  the  application  of  the  nonlinear  WCM,  a greater  number  of  indepen- 
dent tensile  and  compressive  material  properties  are  required  than  for  the  non- 
linear RCM.  Also,  due  to  the  independence  of  tensile  and  compressive  compli- 
ances, restriction  on  the  material's  behavio**  from  the  reciprocal  relations  are 
the  same  for  all  three  energy  functions. 

A short  computer  program,  MULTIW,  was  written  to  investigate  the  weighted 
compliance  model  for  all  three  energy  functions.  MULTIW  is  limited  in  capabil- 
ity to  the  prediction  of  strains  from  an  input  uniform  stress  state.  However, 
MULTIW  has  a degree  of  generality  compatible  with  the  material  models  in  a 
modified  version  of  the  SAAS  III  finite  element  program.  The  stress-strain 
relations  in  MULTIW  can  be  applied  to  an  axi synmetri c body  under  axi symmetric 
load  as  well  as  to  more  general  stress  states  for  isotropic,  transversely  iso- 
tropic, and  orthotropic  materials.  For  the  latter  two  levels  of  anisotropy, 
the  principal  material  (x,  y,  z)  directions  in  Fig.  32  can  exist  at  arbitrary 
orientations  in  the  y-z  plane  to  the  input  stress  {x',  y',  z ' ) directions. 
Because  the  input  stresses  are  constant,  the  principal  stresses  and  their  ori- 
entations do  not  vary  from  one  iteration  to  the  next.  This  program  will  be 
used  in  the  off-axis  strain  response  studies  and  the  biaxial  strain  response 
studies  of  Sections  5 and  6,  respectively. 
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4.5  NONLINEAR  RESTRICTED  COMPLIANCE  MATERIAL  MODEL 


Three  aspects  of  the  nonlinear  restricted  compliance  material  model  (RCM) 
are  discussed:  (1)  the  distinctive  features  of  the  iteration  approach,  (2) 

the  restrictions  imposed  on  the  stress-strain  behavior  that  can  be  treated 
with  the  model,  and  (3)  the  computerized  application  of  the  material  model. 

Only  two  steps  in  the  general  iteration  procedure  described  in  Section  4.2 
are  affected  by  which  multimodulus  procedure  is  used.  These  steps  involve  the 
evaluation  of  the  all-tension  and  all -compress ion  compliance  matrices  in  prin- 
cipal material  coordinates  and  the  evaluation  of  the  multimodulus  compliance 
matrix  in  principal  stress  directions.  The  linear  multimodulus  approach  out- 
lined in  Section  3.3  is  u^ed  to  develop  these  compliance  matrices  for  the  non- 
linear restricted  compliance  material  model.  In  this  approach,  the  symmetry  of 
the  compliance  matrices  about  the  main  diagonal  in  all  coordinate  systems  ic 
ensured  by  the  relation  between  the  material  properties  determined  in  uniaxial 
tension  and  uniaxial  compression.  On  each  iteration  of  the  nonlinear  procedure, 
an  equivalent  linear  (elastic)  system  is  considered.  The  requirement  of  sym- 
metry of  the  compliance  matrix  is  then  related  to  the  existence  of  a potential 
function  for  such  a system. 

The  m_.n  difference  between  the  linear  and  the  nonlinear  RCM  is  the  varia- 
tion of  material  properties  with  strain  energy.  For  the  linear  material  model, 
all  the  input  tensile  and  compressive  material  properties  are  constant.  There- 
fore, the  multimodulus  compliance  matrix  is  affected  only  by  the  signs  of  the 
principal  stresses,  not  their  magnitudes.  In  the  nonlinear  inultimodulus  mate- 
rial model,  the  material  properties  used  to  define  the  all-tension  and  all-com- 
pression .ompliance  matrices  are  functions  of  the  strain  energy.  In  the  first 
iteration,  the  linear  tensile  material  properties  are  used  to  predict  the 
stresses  and  strains.  However,  in  subsequent  iterations,  the  current  value  of 
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the  strain  energy  (or  energies  in  the  case  of  the  divided  energy)  is  used  to 
determine  the  material  properties  in  the  stress-strain  relationship. 


For  the  total  or  weighted  energy  approaches,  the  relationships  between 
material  properties  are  implied  to  be  valid  when  all  the  material  properties 
are  evaluated  at  the  same  strain  energy.  Thus,  restrictions  are  placed  on  the 
uniaxial  stress-strain  behavior  that  can  be  treated  wiLn  the  models.  That  is, 
not  all  materials  have  properties  for  which  the  material  model  restrictions  are 
satisfied. 

Even  more  restrictions  are  imposed  on  the  deformation  behavior  of  the  mate- 
rial being  investigated  if  the  divided  energy  approach  is  used.  In  the  nonlin- 
ear RCM,  the  off-diagonal  compliances  in  corresponding  positions  of  the  all- 
tension  and  all-compression  compliance  matrices  are  required  to  be  identical  at 
the  same  energy  level.  (This  requirement  may  lead  to  limited  applicability  of 
the  model.)  If  the  divided  energy  approach  were  used,  these  off-diagonal  com- 
pliances would  have  to  be  constant  for  all  energy  levels.  A cube  of  an  isotro- 
pic multi  modulus  material  loaded  only  in  two  orthogonal  directions  by  equal  ten- 
sile forces  can  be  used  to  illustrate  the  constancy  requirement  of  the  off-diag- 
onal compliances.  Suppose  two  of  the  three  independent  material  properties  for 
a nonlinear  multi modulus  isotropic  material  were  determined  in  compression  (de- 
spite the  fact  that  we  deal  with  an  all-tension  stress  state).  Then,  the  two 
compression  properties  would  be  used  to  evaluate  all  the  nonzero  off-diagonal 
compliances.  The  compression  component  of  the  total  strain  energy  would  be 
used  to  find  the  current  value  of  the  two  compression  properties.  Because  the 
compression  component  of  total  strain  energy  would  always  be  zero,  these  compli- 
ances are  constant.  Thus,  if  the  material  is  to  satisfy  the  relation 
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at  all  energy  levels,  the  material  property  expressions  for  the  Young's  modulus 
and  Poisson's  ratio  in  tension  could  differ  at  most  in  the  constant  A.  There- 
fore, the  applicability  of  the  nonlinear  RCM  to  actual  deformation  behavior 
must  be  assessed  by  examining  whether  (or  how  well)  the  relations  between  ten- 
sile and  compressive  compliances  are  satisfied. 

A short  computer  program,  MULTIR,  was  written  to  investigate  the  restrict- 
ed compliance  material  model  for  all  three  energy  functions.  This  program  is 
identical  in  form  to  the  MULTIW  computer  program  described  in  Section  4.4  for 
the  weighted  compliance  material  model.  All  conwents  on  the  general  character 
of  the  MULTIW  program  are  also  applicable  to  the  MULTIR  program.  Both  programs 
will  be  used  in  the  off-axis  strain  response  studies  and  the  biaxial  strain 
response  studies  of  Sections  5 and  6,  respectively. 

4.6  SUMMARY 

The  nonlinear  material  model  in  Section  2 is  combined  with  each  of  the 
linear  multimodulus  material  models  in  Section  3 to  obtain  the  two  nonlinear 
multimodulus  material  models  which  are  described  in  this  Section.  The  use  of 
three  different  energy  functions  is  investigated  for  both  of  the  two  models. 
Restrictions  on  the  stress-strain  behavior  that  can  be  treated  are  discussed. 

The  actual  stress-strain  behavior  of  a nonlinear  multimodulus  material 
would  probably  have  the  general  form  of  the  dashed  curve  in  Fig.  5.  In  the 
linear  multimodulus  models,  a bilinear  approximation  of  the  actual  stress- 
strain  behavior  is  used.  This  bilinear  approximation  is  represented  by  the 
solid  lines  in  Fig.  5.  The  discontinuity  in  slope  at  the  origin  for  the  bi- 
linear approximation  would  probably  not  occur  for  an  actual  material.  Instead, 
a nonlinear  transition  rec:on  would  be  expected.  This  transition  region  can 
be  represented  with  the  nonlinear  multimodulus  material  models  by  use  of  the 
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same  value  of  the  material  property  constant  A,  but  different  values  of  B and 
C,  in  the  tension  and  compression  equations  for  a material  property.  For 
problems  in  which  the  stress-strain  behavior  in  the  transition  region  is  not 
important,  a better  representation  of  the  actual  material  property -strain 
energy  variations  might  be  achieved  with  different  values  of  A in  the  tension 
and  compression  equations  for  a material  property. 

The  applicability  and  accuracy  of  the  nonlinear  multimodulus  material 
models  are  investigated  in  the  off-axis  strain  response  and  biaxial  strain 
response  studies  of  Sections  5 and  6,  respectively. 


5.  COMPARISON  OF  PREDICTED  AND  MEASURED  STRAIN  RESPONSE 


UNDER  UNIAXIAL  OFF-AXIS  LOADING 

5.1  INTRODUCTION 

The  two  nonlinear  multimodulus  material  models  described  in  Section  4 are 
used  to  predict  the  uniaxial  deformation  behavior  of  ATJ-S  graphite  when  the 
load  is  applied  at  an  angle  a to  the  principal  material  directions  (see  Fig. 32). 
These  predictions  are  compared  to  the  strain  response  measured  by  Jortner  [30], 
The  comparisons  between  predicted  and  measured  strains  are  used  to  evaluate  how 
well  the  uniaxial  deformation  behavior  can  be  predicted  when  (1)  the  load  is 
applied  in  the  principal  material  directions  and  (2)  the  load  is  applied  in 
other  than  principal  material  directions.  In  addition,  the  application  of  the 
models  and,  to  some  extent,  the  contrast  between  them  is  demonstrated. 

Before  presentation  and  discussion  of  the  predicted  and  measured  strain 
results,  the  pertinent  features  of  the  test  specimens  and  the  test  procedures 
are  described. 

5.2  DESCRIPTION  OF  TEST  SPECIMENS  AND  PROCEDURES 

The  off-axis  specimens  were  obtained  from  ATJ-S  graphite  billet  16K9-27  by 
the  scheme  shown  in  Figs.  33-35.  The  nominal  directions  of  the  planes  of  iso- 
tropy of  the  transversely  isotropic  material  are  identified  with  light  lines. 

The  billet  cylindrical  coordinate  system,  shown  in  Fig.  33,  is  nominally  aligned 
with  the  principal  material  directions.  That  is,  the  principal  material  direc- 
tions are  not  necessarily  the  same  throughout  the  billet  because  of  nonuniform- 
ities,  e.g.,  ncnuniform  compaction  or  nonuniform  particle-binder  mixture,  in 
the  billet  manufacturing  process.  The  6-in.  long  tension  test  specimens  -..ith 
a 1/4-in.  square  gage  section  were  cut  from  two  of  the  slabs  from  the  billet 
as  schematically  shown  in  Fig.  34.  The  2-in.  long  compression  test  specimens. 
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FIGURE  33  RELATION  OF  TEST  SPECIMEN  SLABS  TO  BILLET 


which  were  of  the  characteristic  dogbone  shape,  nave  a 1/2-inch  diameter  cylin- 
drical gage  section.  These  compression  specimens  were  cut  from  a slab  of  the 
billet  as  schematically  shown  in  Fig.  35.  At  least  three  tension  and  compres- 
sion test  specimens  were  cut  from  the  slabs  at  each  of  four  angular  orienta- 
tions, 0°,  45°,  70°,  and  90°,  from  principal  material  directions. 

The  positions  of  the  axial  and  transverse  strain  gages  on  the  gage  section 
of  the  tension  and  compression  specimens  are  indicated  in  Figs.  36  and  37.  The 
figures  are  not  adequate  for  precise  definition  of  the  instrumentation  of  each 
specimen.  (A  detailed  description  of  the  strain  gages  on  each  specimen  is  given 
in  Tables  Xill  and  XIV  of  P.ef.  /.)  However,  sufficient  information  is  in  Figs. 
36  and  37  for  consistent  comparison  of  predicted  and  measured  strains. 

The  (x'-y'-z1)  coordinate  system  used  in  the  presentation  of  strain  results 
is  shown  in  Figs.  36  and  37.  Note  that  the  x-y-z  coordinate  system  is  aligned 
with  the  nominal  principal  material  directions  of  the  material  and  that  the  x 
and  y coordinates  are  in  the  plane  of  isotropy  of  the  material.  The  angle  i is 
the  angle  in  the  y-z  plane  between  the  load  axis  (z’)  and  the  across  grain  (z) 
principal  material  direction.  For  example,  the  specimen  is  loaded  uniaxially 
in  the  across  grain  (z)  direct.-' on  when  j = 0C  and  in  the  with  grain  (y ) direc- 
tion when  j.  = 90°.  Strains  were  always  measured  in  the  x'-y’-z1  coordinate 
directions. 

Acoustic  velocity  measurements  and  shear  coupling  effects  on  tne  specimens 
are  reported  by  Jortner  [7],  Acoustic  velocity  measurements  we’-e  made  tc  deter- 
mine variations  between  local  and  nominal  principal  material  directions  due  to 
billet  nonuni formi ties . As  a result  of  this  investigation,  ,’ortner  speculates 
that,  for  the  45°  tension  specimens  and  the  70°  compression  specimens,  the 
angles  between  principal  material  and  load  directions  may  be  Closer  to  40c  and 
80°,  respectively.  Analytical  studies  of  shear  coupling  effects  on  the  stress 
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FIGURE  36  TENSION  BAR  GAGE  SECTION 
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FIGURE  37  COMPRESSION  ROD  GAGE  SECTION 


distributions  in  the  off-axis  specimens  are  also  presented  in  Ref.  7.  For 
both  the  tension  and  compression  specimens,  the  effect  of  sne?:  coupling  on 
the  stress  distribution  across  the  width  of  the  gage  sections  is  estimated 
to  be  sma 1 1 . 

The  specimens  were  loaded  on  ar.  Instron  test  machine.  The  load  was  intro- 
duced in  the  tension  specimens  with  doublers  bonded  to  the  specimens.  ;'hese 
doublers  were  loaded  tnrough  pin  attachments.  The  compression  specimens  were 
loaded  in  a test  fixture  with  an  estimated  friction- induced  error  in  the  re- 
corded loaa  of  less  than  2 [ 3 j j . 

5.3  COMPARISONS  OF  PREDICTED  AND  MEASURED  STRAIN  RESPONSE 

Discussion  of  correlation  between  predicted  and  measured  strain  response 
is  dividec  into  five  sections.  In  Section  5.3.1,  the  data  base  material  prop- 
erties are  discussed.  The  general  method  used  to  present  and  compare  the  pre- 
dicted and  measured  strain  response  under  tensile  and  com;  res  si ve  leads  is 
described  in  Section  5.3.2.  The  tensile  strain  correlations,  «h''"h  are  the 
same  for  both  the  nonlinear  ViCM  and  RCM,  are  presented  m Secfc'-  r .3.3.  The 
compressive  strain  correlations  are  discussed  for  the  nonlinear  ,r  Section 
5.3.4.  Finally,  the  compressive  strain  correlations  f_>r  ih-e  nonlinea-  RCM  are 
described  in  Section  5.J.5. 

5.3.1  Materia i Properties 

The  first  step  in  application  of  the  new  material  model:  in  strain  corre- 
lation studies  is  to  develop  a set  of  material  property  , tan is  f^r  the  mate- 
rial property-strain  energy  equations  from  appropriate  uniaxial  te'-t  results. 

The  material  property  constants  used  in  the  present  strain  correct  ions  .ere 
determined  directly  from  ’ortner's  experimental  results  [30j.  Ihus,  extraneous 
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variations  between  predicted  and  measured  strains  are  not  introduced  because 
of  bi llet-to-bi 1 let  variations  in  material  properties  or  differences  in  test 
procedures.  The  material  property  constants  in  the  billet  coordinate  system 
(principal  material  coordinates)  are  given  in  Table  8. 

Different  independent  material  properti. i are  required  for  each  of 
the  two  nonlinear  multi  modulus  material  models.  For  a transversely  isotropic 
multimoduius  material,  six  independent  material  properties  are  needed  to  apply 
the  restricted  compliance  material  model.  The  predictions  of  the  RCM  are  de- 
pendent on  which  material  properties  are  selected  for  the  uniaxial  baseline 
values.  Of  course,  any  set  of  six  material  properties  could  be  used  and  the 
same  predictions  would  result  if  the  restrictions  of  the  material  model  were 
perfectly  satisfied  by  the  stress-strain  behavior  of  the  material.  For  trans- 
versely isotropic  multimoduius  materials,  ten  independent  material  properties' 
are  required  to  apply  the  weighted  compliance  material  model.  The  input  mate- 
rial property  sets  investigated  are: 

Weighted  Compliance  Material  Model 

Ert’  Ezt’  vrt-t*  '2et*  Erzt’ 

Erc’  Ezc’  vr9c’  'zvc'  E rzc 
Restricted  Compliance  Material  Model 

Property  Set  E^:  Ert,  >r9f  “zef  Em'  Ezc 


Property  Set  E 


F f v V £'  F 
rt*  ~zt’  VrSt’  ret’  trzt’  rc 


5.3.2  Method  of  Presentation  of  Results 

Strain  predictions  with  tne  new  material  models  are  shown  in  Figs.  38- 
53  as  symbols  next  to  the  curves  for  experimental  results.  The  experimental 
results  were  obtained  from  Jortner's  Figs.  4 thru  7 and  9 thru  12  [30].  The 
predicted  strains  were  obtained  with  the  MULTIR  and  MULTIW  computer  programs. 
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TABLE  8 CONSTANTS  IN  MATERIAL  PROPERTY  EQUATIONS  FOR 
16K9  BILLET  DATA  BASE*  AT  70° F 


MATERIAL  PROPERTY 


1 .37xlO°psi 

.1 

.1607 

1. 385x1 06psi 
1.5152xl06psi 
1.1 538xl06ps i 
.09 
.1 

1 .2626xl06psi 


B 

C 

.101576 

.445843 

.231769 

.255095 

.376247 

.169649 

.445588 

.157836 

.132548 

.365323 

.0609584 

.646868 

.140150 

.370244 

0 

1 

0 

1 

.122519 

.394504 

*Experimenta1  Data  Due  to  Jortner  [30],  Figs.  4,  5,  7,  9,  10,  and  12 


Because  a portion  of  the  experimental  results  was  used  to  develop  the  mate- 
rial property  constants  input  to  these  computer  programs,  different  categories 
of  strain  predictions  must  be  discussed.  The  material  property  constants 
were  determined  by  a simultaneous  solution  technique  from  data  at  three  points 
on  the  experimental  curves.  Predictions  of  strains  for  these  data  base  points 
are  indicated  with  dots  (fill ed--* n circles),  whereas  other  predictions  of  the 
material  model  are  indicated  with  open  circles.  All  the  predictions  along 
stress-strain  curves  used  to  define  material  property-strain  energy  relation- 
ships are  distinguished  from  other  predictions  by  the  appearance  of  the  mate- 
rial property  designation  next  to  its  defining  curve.  For  example,  the  experi-  ; 

f 

mental  results  in  Fig.  38  were  used  to  define  the  constants  in  the  equation 
relating  material  properties  E^  and  v to  the  strain  energy.  Predictions 

for  the  actual  points  used  in  the  derivation  of  the  material  property  constants  j 

« 

are  indicated  with  dots. 

In  the  discussion  of  results,  approximate  percentage  differences  between 
predicted  and  measured  strains  are  computed  from 
Predicted  - Measured 

As  = x 100%  I 

Measured  (31) 

When  more  than  one  experimental  curve  is  presented  for  stress-strain  behavior 
other  than  data  base  behavior  (curves  used  to  compute  the  material  property 
constants) , the  percentage  differences  are  calculated  between  the  predicted  value 
and  the  closest  experimental  results  (a  minor  bius).  For  example,  the  perceutace  | 

difference  for  the  transverse  strain  prediction  in  Fig.  39  at  4000  psi  would  be 
computed  with  the  strain  from  the  left-most  experimental  curve.  When  the  predic- 
tion falls  between  experimental  curves,  the  percentage  differences  are  regarded 
as  zero. 
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5.3.3  Tensile  Strain  Correlations 

The  predictions  with  both  material  models  for  each  of  the  three  energy 
approaches  are  identical  for  the  tension  loading  cases  in  Figs.  38-41.  The 
predictions  are  identical  because  of  the  formulations  of  the  two  material  models 
and  the  use  of  five  tension  properties  in  both  property  sets  for  the  restricted 
compliance  material  model.  The  curves  used  to  define  the  material  property- 
strain  energy  relationships  for  the  five  tension  properties  common  to  all  three 
property  sets  (WCM  and  both  RCM)  are  shown  in  Figs.  38,  39,  and  41.  The  largest 
variation  between  predicted  and  measured  strains  which  occurs  for  a curve  used 
to  establish  a direct  moduli-strain  energy  relationship  is  5%.  The  predictions 
for  curves  in  Figs.  38  and  41  which  were  used  to  determine  material  property 
constants  for  \>z^  and  v , respectively,  show  equally  good  (within  about  5%) 
correlation  except  for  the  1000  psi  and  2000  psi  stress  level  predictions  in 
Fig.  41.  The  strains  predicted  for  1000  and  2000  psi  stress  levels  differ 
from  the  measured  strains  for  the  v curve  by  -19%  and  -14%,  respectively. 
These  variations  are  not  considered  significant  because  they  occur  at  low  strain 
levels  which  are  generally  not  critical  in  structural  design  analysis.  Varia- 
tions between  predicted  and  measured  strains  for  1000  psi  and  2000  psi  stress 
states  on  the  vrgt  curve  might  be  attributed  to  a poor  approximation  of  the 
material  property-strain  energy  relationships  at  low  energy  levels.  A low-fit 
approximation  as  shown  in  Fig.  22  might  result  in  improved  strain  correlation 
for  these  low  stress  levels.  A visual  inspection  of  the  strain  correlations 
for  the  vrQt  curve  in  Fig.  41  reveals  that  the  variations  at  1000  psi  and  2000 
psi  are  not  of  practical  importance  because  the  strains  themselves  are  quite 
small.  The  use  of  | REL  AUj  = .5%  as  a convergence  test  value  in  the  MULTIR 
and  MULTIW  computer  programs  is  supported  by  the  excellent  strain  correlations 
found  for  the  input  data  points  (filled-in  circles)  shewn  in  Figs.  38,  39,  and 
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41.  Good  strain  correlations  are  found  for  the  predictions  {open  circles)  of 
points  on  the  curves  used  to  evaluate  the  material  property  constants  for  the 
five  tension  properties.  Thus,  the  general  material  property  equation,  Eq.  (5) 
can  be  used  to  represent  the  material  property-strain  energy  variations  for  the 
tension  material  properties  of  ATJ-S  graphite  at  70°F.  The  simultaneous  solu- 
tion procedure  used  to  evaluate  the  material  property  constants  B and  C is 
also  generally  substantiated  by  the  comparisons. 

Four  sets  of  strain  correlations  for  the  tensile  loads  are  not  directly 
associated  with  input  material  properties:  (1)  transverse  strain  zy,  at  a=45° 

(cig.  39),  (2)  axial  strain  e2i  at  a=70°  (Fig.  40),  (3)  transverse  strain  ey> 
at  of 70°  (Fig.  40),  and  (*)  transverse  strain  cyt  at  a=90°  (Fig.  41).  Good 
strain  correlation  was  found  for  £z.  at  u=70°  and  eyi  at  cf90°,  whereas  poor 
agreement  between  predicted  and  measured  strains  exists  for  cy,  at  <x=45°  and 
70°.  Specifically,  at  3000  psi,  the  predicted  axial  stress  at  a=70°  differs 
from  the  measured  results  by  approximately  2%.  The  estimated  percentage  dif- 
ference for  the  transverse  strain  cy,  at  a=90e  is  6%.  The  good  agreement  be- 
tween the  predicted  and  measured  transverse  strains  eyt  in  Fig.  41  is  an  indi- 
cation that  the  linear  reciprocal  relations  between  compliances  are  applicable 
to  the  nonlinear  case.  At  a stress  of  3000  psi,  transverse  strain  predictions 
at  of45°  {Fig.  39)  and  a=70°  (Fig.  40)  differ  from  the  measured  results  by 
approximately  54%  and  22%,  respectively.  No  suitable  explanation  exists  for 
the  disparity  between  predicted  and  test  values  of  the  transverse  strain  at 
a=45°  and  70°.  However,  an  anomalous  feature  of  the  cf45°  specimens,  namely, 
a variation  in  orientation  of  the  principal  material  directions,  was  uncovered 
in  the  acoustic  velocity  investigation.  Also  worthy  of  note  is  the  fact  that 
if  the  strain  predictions  for  exi  are  shown  on  Figs.  39  and  40,  then  the  mea- 
sured transverse  strains  are  between  the  exi  and  ev»  predictions. 


5.3.4  Compressive  Strain  Correlations  for  Nonlinear  Weighted  Compliance  Model 
The  measured  strains  for  specimens  under  uniaxial  compression  loads  are 
shown  with  predictions  from  the  nonlinear  weighted  compliance  model  in  Figs. 
42-45.  For  these  loading  conditions,  the  WCM  predictions  are  identical  for  all 
energy  approaches.  In  all  instances,  good  correlation  between  predicted  and 
measured  strains  is  found.  The  largest  percentage  deviation  at  compressive 
stress  levels  with  magnitudes  greater  than  3000  psi  is  12%  and  occurs  for  the 
transverse  strain  at  a=70°  (Fig.  44)  at  a stress  of  -5000  psi.  This  transverse 
strain  correlation  would  be  improved  if  a were  regarded  as  80°  as  found  in  the 
acoustic  velocity  investigation.  The  curves  used  to  define  the  material  proper- 
ty-strain energy  relationships  for  the  five  independent  compressive  properties 
in  the  WCM  data  set  are  shown  in  Figs.  42,  43  and  45.  The  largest  variation 
between  predicted  and  measured  strains  which  occurs  for  a curve  used  to  find 
the  material  property  constants  is  8%.  Thus,  the  general  material  property 
equation,  Eq.  (5),  can  be  used  to  represent  the  material  property-strain  energy 
variations  for  the  compression  material  properties  of  ATJ-S  graphite  at  70°F. 

The  procedure  used  to  evaluate  the  material  property  constants  in  the  general 
equation  is  also  substantiated  by  the  results.  Some  of  the  deviation  between 
the  experimental  axial  strain  curves  in  Fig.  42  (also  Figs.  46  and  50)  might 
be  due  to  a calibration  error  for  the  lowest  curve  (Jortner  [30]  indicated  a 
possible  calibration  error  with  a comment  by  the  curve). 


5.3.5  Compressive  Strain  Correlations  for  Nonlinear  Restricted  Compliance  Model 
The  restricted  compliance  material  model  is  applied  for  compressively 
loaded  specimens  to  obtain  the  predicted  results  in  Figs,  46-53.  Results  for 
material  property  setr  Ezc  and  Erc  are  displayed  in  Figs.  46-49  and  Figs.  50-53, 
respectively.  Recall  that  the  two  material  property  sets  differ  by  only  one 
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independent  material  property.  The  material  property  E2C  used  in  the  first 
property  set  is  replaced  by  Erc  in  the  second  property  set.  All  three  energy 
approaches  are  investigated  for  material  property  set  Ezc.  Only  the  total  en- 
ergy approach  is  examined  for  the  econd  property  set  Erc.  The  dependence  of 
the  predictions  cf  the  RCM  on  the  selection  of  the  baseline  material  properties 
is  demonstrated  by  the  results.  For  example,  at  5000  psi,  the  axial  strain 
prediction  for  the  total  energy  case  in  Fig.  46  is  12%  smaller  than  the  same 
prediction  in  Fig.  50.  Similar  deviations  (8-10%)  in  axial  strains  occur  when 
the  predictions  in  Figs.  51,  52,  and  53  are  compared  with  the  predictions  in 
Figs.  47,  48,  and  49,  respectively.  The  differences  in  transverse  strain  pre- 
dictions .for  the  total  energy  approach  with  the  two  different  property  sets  are 
small  (less  than  5%)  because  the  off-diagonal  compliances  are  determined  from 
the  same  independent  material  properties  (tension  properties)  in  both  cases. 
Differences  in  transverse  strains  predictions  occur  because  the  material  prop- 
erties are  determined  at  differen.  .1  energy  levels.  Because  the  corre- 
sponding axial  strain  predictions  wiui  the  different  data  sets  are  not  the 
same,  the  relations  between  tension  and  compression  compliances  in  the  RCM 
may  not  be  strictly  applicable  to  ATJ-S  graphite. 

Although  the  axial  strains  in  Figs.  46-49  agree  relatively  well  (within 
9%)  with  test  results  for  the  total  and  weighted  energy  approaches,  the  trans- 
verse strain  correlation  is  generally  poor.  The  largest  deviation  in  axial 
strains,-9%,  occurs  for  the  a=90°  case  (Fig.  49).  At-5000  psi,  all  transverse 
strain  predictions  for  the  total  and  energy  approaches  in  Figs.  46-49  differ 
from  measured  results  by  44%  to  55%.  The  relationships  between  diagonal  com- 
pliances (direct  moduli)  for  the  all-tension  and  all-compression  stress  states 
in  the  nonlinear  RCM  might  be  applicable  to  ATJ-S  graphite,  but  those  between 
off-diagonal  Compliances  (in  which  both  direct  moduli  and  Poisson's  ratios 
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appear)  are  not  applicable.  The  restrictions  on  the  off-diagonal  compliances 
in  the  nonlinear  ROM  may  be  especially  stringent  for  nonlinear  stress-strain 
behavior  of  the  material.  When  all  material  properties  are  evaluated  at  the 
same  energy  level,  the  compliance  relationships  should  be  valid.  In  the  70°F 
experimental  results  in  Refs.  30  and  14,  the  tension  Poisson's  ratios  gener- 
ally decrease  with  increasing  energy  level,  whereas  Poisson's  ratios  in  com- 
pression usually  remain  almost  constant  or  increase.  Thus,  the  variations  for 
Poisson's  ratios  implied  in  the  nonlinear  RCM  may  not  be  representative  of  the 
actual  variations.  When  the  divided  energy  approach  is  used  in  the  nonlinear 
RCM  instead  of  the  total  or  weighted  energy  approaches,  the  transverse  strain 
correlation  is  somewhat  improved  (the  largest  deviation  is  about  28%),  but  the 
axial  strain  correlation  is  much  poorer  (the  largest  deviation  is  approximately 
35%).  With  the  divided  energy  approach,  the  implied  off-diagonal  compliances 
are  constant  on  successive  iterations  for  the  loading  and  independent  material 
properties  investigated.  The  improved  correlation  for  the  transverse  strain 
results  from  the  larger  strain  predictions  in  the  divided  energy  approach  than 
in  the  total  or  weighted  energy  approaches.  These  predictions  occur  because 
the  implied  off-diagonal  compliances  remain  constant  instead  of  decreasing  with 
increasing  energy  level.  The  implied  diagonal  compliances  in  the  divided  energy 
approach  seem  to  be  inappropriate  because  of  the  poor  axial  strain  correlations. 

5.4  SUMMARY 

The  strain  predictions  with  both  the  nonlinear  restricted  compliance  and 
weighted  compliance  material  models  are  compared  with  the  experimental  strain 
response  obtained  from  uniaxial  loading  of  specimens  with  principal  material 
directions  at  foir  different  orientations  with  respect  to  the  load  axis.  The 
best  overall  correlation  between  predicted  and  measured  strains  occurs  for  the 
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nonlinear  weighted  compliance  material  model.  Note  that  the  nonlinear  weighted 
compliance  material  model  results  are  the  same  for  all  energy  approaches  for 
these  uniaxial  stress  conditions.  That  is,  because  the  stresses  are  either  all 
tension  or  all  compression,  the  three  energies  used  in  each  of  the  three  energy 
approaches  are  identical.  Moreover,  the  compliances  are  all  determined  at  the 
same  energy  level.  In  contrast,  for  the  nonlinear  restricted  compliance  model, 
the  total  and  weighted  energies  are  the  same  but  the  divided  energy  approach 
leads  to  an  undesirable  consequence.  Specifically,  in  an  all-compression  stress 
state,  the  compression  energy  is  the  total  energy  and  the  tension  energy  is 
zero;  accordingly,  the  compliances  are  determined  from  a mixed  set  of  nonlinear 
compression  properties  and  linear  (elastic)  tension  properties.  Thus,  differ- 
ences between  measured  strains  and  strains  predicted  from  the  nonlinear  restricted 
compliance  material  model  are  attributed  mainly  to  the  relations  between  tension 
and  compression  compliances  used  in  the  model.  The  stress-strain  behavior  of 
ATJ-S  graphite  evidently  does  not  closely  satisfy  the  restrictions  between  com- 
pliances in  the  nonlinear  RCM. 

For  all  direct  comparisons  between  predicted  and  measured  strains  for  the 
data  base  material  properties,  excellent  correlation  is  obtained.  Thus,  the 
general  material  property  equation,  Eq.  (5),  can  be  used  to  represent  the  mate- 
rial property-strain  energy  variation  for  both  tension  and  compression  proper- 
ties of  ATJ-S  graphite  at  70°F.  The  procedure  used  to  evaluate  the  material 
property  constants  B and  C is  also  substantiated  by  the  comparisons. 

The  applicability  of  the  nonlinear  multi  modulus  material  models  to  the 
prediction  of  the  biaxial  strain  response  of  ATO-S  graphite  is  investigated  in 
Section  6. 
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6.  COMPARISON  OF  PREDICTED  AND  MEASURED  STRAIN  RESPONSE 


UNDER  BIAXIAL  LOADING 


6.1  INTRODUCTION 

The  two  nonlinear  multi  modulus  material  models  discussed  in  Section  4 are 
used  to  predict  the  biaxial  strain  response  of  ATJ-S  graphite  in  this  section. 
Strain  predictions  of  both  material  models  are  compared  with  strains  measured 
by  Jortner  [6-9,  30]  in  room  temperature  tests.  Strain  predictions  of  the  non- 
linear material  model  of  Section  2 are  also  compared  with  the  biaxial  tension 
strain  response  measured  at  room  temperature.  Thus,  an  assessment  of  the  ba- 
sic nonlinear  approach  without  the  complications  of  the  multimodulus  proce- 
dures can  be  made.  Biaxial  fracture  investigations  reported  by  Jortner  [6-9] 
are  another  source  of  experimental  data  which  can  be  used  to  assess  the  appli- 
cability of  the  nonlinear  multimodulus  material  models  to  prediction  of  ATJ-S 
graphite  stress-strain  behavior.  Strains  at  failure  measured  in  at  room  tem- 
perature tests  are  compared  with  predicted  strains.  The  correlation  between 
the  predicted  strain  response  with  the  nonlinear  restricted  compliance  mate- 
rial model  and  the  experimental  strain  response  at  2000°F  determined  by  Jort- 
ner [7,  8]  is  also  examined. 

The  predicted  and  measured  strain  results  for  the  room  temperature  bi- 
axial strain  study,  the  room  temperature  fracture  strain  study,  and  the  2000°F 
biaxial  strain  study  are  presented  in  Sections  6.2,  6.3,  and  6.4,  respectively. 
The  pertinent  features  of  the  test  specimens,  test  procedures,  and  data  base 
material  properties  are  discussed  in  each  section.  Effects  in  the  strain  cor- 
relations from  differences  in  specimen  origin  (different  billets),  friction  on 
the  measured  axial  loads,  and  circumferential  stress  gradients  in  the  biaxial 
specimens  are  examined.  Adjustments  made  for  these  factors  are  briefly  de- 
scribed. 
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In  the  discussion  of  results  In  this  section,  percentage  deviations  be- 
tween predicted  and  measured  strains  are  frequently  used.  Recall  from  the  dis- 
cussion of  Section  2 that  percentage  deviation  is  defined  as  the  distance  be- 
tween the  experimental  (baseline)  and  predicted  strains  on  the  biaxial  strain 
plot  divided  by  the  length  of  a ray  from  the  origin  to  the  experimental  strain 
values.  The  percentage  deviation  is  schematically  depicted  in  Fig.  23  and  is 
calculated  from 


% Deviation  = 


\ + 


>2  "1 1/2 


x 100% 


where  the  subscripts  have  the  following  connotations: 

a - axial  direction 
c - circumferential  direction 
m - measured  values 
p - predicted  values 

The  percentage  deviation  is  often  preceded  by  a + or  - sign  to  denote  whether 
the  predicted  strains  are  an  overestimate  (+)  or  an  underestimate  (-)  of  the 
measured  strains.  The  measured  strains  are  only  approximate  since  most  of 
the  experimental  results  discussed  in  this  section  were  obtained  from  graphi- 
cal results  instead  of  tabulated  results.  Therefore,  the  percentage  devia- 
tion calculations  are  approximate. 


6.2  ROOM  TEMPERATURE  STRAIN  CORRELATIONS 

Jortner's  measured  biaxial  strain  results  at  room  temperature  are  com- 
pared with  nonlinear  WCM  and  nonliner  RCM  predictions  in  this  section.  Bi- 
axial tension  strain  response  is  also  compared  with  nonlinear  material  model 
predictions  to  obtain  an  assessment  of  the  nonlinear  approach  without  the  com- 
plications of  the  multimodulus  models. 
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A typical  hollow  biaxial  test  specimen  is  shown  in  Fig.  54  (from  Fig.  9, 
Ref.  7).  Internal  pressure  is  applied  to  the  specimen  by  a pressurizing  fluid 
contained  in  a thin  rubber  bladder.  Axial  loads  (either  tension  or  compression) 
are  applied  to  the  specimen  at  the  lips  on  the  specimen  ends  through  a fixture 
designed  to  restrict  bending.  Thus,  tests  can  be  performed  in  two  quadrants  of 
biaxial  loading.  Specimens  were  loaded  to  several  different  ratios  of  axial 
stress  to  circumferential  stress  (o  :o  ).  These  nominal  stress  ratios  were 

d C 

approximately  constant  during  loading  to  a maximum  principal  thin-wall  stress 
of  3550  psi.  Strains  were  measured  with  strain  gages  on  the  outer  surface  at 
the  midlength  of  the  specimens.  The  nominal  principal  material  directions  are 
aligned  with  the  load  axes  for  all  specimens.  The  billet  coordinate  system  is 
appropriate  to  use  to  describe  material  properties  for  the  biaxial  test  speci- 
mens since  the  billet  axes  are  also  aligned  with  the  nominal  principal  material 
directions.  The  specimen  axis  (axial  load  direction)  is  aligned  with  the  billet 
axis  (z  - direction). 


CIRCUMFERENTIAL 
4.00 — 


1.01  1.68 


DIMENSIONS  IN  INCHES 


FIGURE  54  HOLLOW  GRAPHITE  BIAXIAL  TEST  SPECIMEN  WITH  STRAIN  GAGES 
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6.2.1  Strain  Correlations  for  Nonlinear  Material  Model  Predictions 

Strain  correlations  are  performed  to  determine  the  applicability  of  the 
nonlinear  material  model  in  Section  2 to  prediction  of  biaxial  stress-strain 
behavior  ATJ-S  graphite.  The  nonlinear  multimodulus  material  models  in 
Section  4 effectively  degenerate  to  the  nonlinear  material  model  either  if 
the  material  property  constants  are  identical  for  each  independent  tension 
and  compression  property  and  the  total  energy  approach  is  used  or  if  an  all- 
compression  or  an  all-tension  stress  state  exists  (with  some  obvious  data  input 
restrictions  for  the  RCM  model).  In  this  section,  biaxial  tension  stress  states 


are  examined  which  result  from  internal  pressure  and  axial  tension  on  the  speci- 


men in  Fig.  54.  Actually,  the  induced  stresses  are  not  all  tension,  but  the 

i 

radial  compressive  stress  is  very  small  in  comparison  to  the  circumferential 
and  axial  tensile  stresses.  Thus,  the  applicability  of  the  nonlinear  material1 
model  can  be  determined  without  the  complications  of  the  multimodulus  features 
of  the  WCM  and  RCM  models. 

To  evaluate  the  effects  of  stress  gradients  in  the  specimens  on  the  theo- 
retical-experimental comparisons,  the  grid  with  355  finite  elements  (provided 
by  Jortner)  in  Fig.  55  was  used  to  model  each  of  the  specimens.  Only  half  of 
an  actual  test  specimen  need  be  modeled  because  of  the  symmetry  of  the  speci- 
mens and  applied  loads.  This  grid  was  used  in  a modified  version  of  the  SAAS 
III  finite  element  computer  program  with  the  new  nonlinear  material  model.  The 
predicted  stresses  and  strains  are  fairly  uniform  throughout  the  thinned  down 
portion  of  the  specimen.  However,  there  is  a detectable  stress  gradient  throuyh 
the  specimen  wall  due  to  bending.  Thus,  the  specimen  in  Fig.  54  is  a reasonably 
good  engineering  test  specimen. 

The  measured  strains  and  the  strains  predicted  with  the  new  nonlinear  ! 
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material  model  in  the  SAAS  III  program  are  shown  in  Fig.  56.  The  measured  | 

strain  values  for  a maximum  principal  thin-wall  stress  of  3550  psi  are  Jort-  1 
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FIGURE  55  FINITE  ELEMENT  MOOEL  OF  HOLLOW  GRAPHITE  BIAXIAL  TEST  SPECIMEN 


ner's  results  £8,  30.1  on  the  outside  of  the  specimen  midlength.  In  all  but  one 
case,  the  experimental  results  shown  in  Fig.  56  are  for  specimens  from  a single 
billet  (ICO— 15) . The  measured  strains  for  the  equal  biaxial  tension  stress  state 
are  from  a billet  3R9-33  specimen.  For  the  nominal  stress  ratios  (1:1.26)  and 
(-1:1),  two  data  points  with  slightly  different  actual  stress  ratios  are  pre- 
sented. 

Strain  predictions  for  the  nominal  stress  ratios  at  a maximum  principal 
th in-wall  stress  of  3550  psi  are  shown  in  Fig.  56  for  the  nonlinear  material 
model.  For  the  convergence  criterion  used  (relative  change  in  total  strain 
energy  .005),  the  procedure  converged  in  either  five  or  six  Iterations. 

When  the  finite  element  procedure  is  used,  stress  and  strain  gradients  can 
be  estimated  in  the  specimens.  The  strain  predictions  shown  are  extrapolated 
to  the  outer  surface  of  the  specimens  (where  the  strains  were  measured)  from 
the  center  of  the  element  just  inside  the  outer  surface. 

The  SAAS  II  predictions  in  Ref.  8 for  a maximum  principal  thin-wail 
stress  of  35G0  psi  are  repeated  in  Fig.  56  to  allow  comparison  of  the  new 
nonlinear  material  model  with  the  material  model  now  used  in  reentry  vehicle 
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nosetip  stress  analysis  (SAAS  II  and  SAAS  III),  In  the  material  model  of 
SAAS  II  [32],  a bilinear  representation  of  the  uniaxial  stress-strain  data  is 
used  in  conjunction  with  a "normalized  effective  stress"  function  to  describe 
nonlinear  orthotropic  deformation  behavior.  There  is  no  explicit  assumption 
of  a zero  plastic  volume  change  in  SAAS  II.  Instead,  the  stress  function  that 
is  used  in  SAAS  II  and  the  unaltered  version  of  SAAS  III  [24]  reduces  for  iso- 
tropic materials  to  the  von  Hises  distortional  energy  relation.  However,  that 
stress  function  is  an  arbitrary  expression  that  cannot  be  derived  (that  is,  it 
has  no  basis  on  theoretical  grounds).  Despite  that  necessary  condition  of  re- 
duction to  a known  relation  (which  is  appropriate  for  isotropic  metals  but 
has  not  been  shown  to  apply  to  granular  materials),  the  sufficiency  argument 
for  orthotropic  materials  has  not  been  satisfied.  Thus,  the  applicability  of 
the  nonlinear  material  models  in  SAAS  II  and  the  SAAS  III  computer  programs 
to  prediction  of  stress-strain  behavior  for  ATJ-S  graphite  is  questionable. 

In  addition,  the  multimodulus  character  of  the  material  cannot  be  represented 
with  SAAS  II  material  model.  The  SAAS  III  program  does  have  the  linear  WCM 
discussed  in  Section  3.2.  However,  the  nonlinear  and  multimodulus  material 
models  cannot  be  used  simultaneously  in  SAAS  III.  Accordingly,  any  agreement 
obtained  between  measured  strains  and  strains  predicted  with  the  SAAS  II  (or 
SAAS  III)  material  model  is  both  coincidental  and  fortuitous. 

To  reduce  effects  caused  by  differences  in  billet  stiffness,  specimen 
design,  and  test  procedures  in  the  evaluation  of  the  nonlinear  material  model, 
uniaxial  tests  of  the  tubular  specimens  from  billet  1C0-15  were  used,  where 
possible,  to  establish  the  material  property  constants.  The  material  prop- 
erty constants  [A,  B,  C and  U0  of  Eq.  (5)]  for  the  five  independent  material 
properties  are  given  in  billet  coordinates  in  Table  9.  The  following  five 
independent  tensile  material  properties  are  used  as  a data  base  for  the  non- 
linear material  model:  Ert,  Ezt,  v^,  vzQt,  and  Gr2t*  Differences  between 
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TABLE  9 

ICO  DATA  BASE  AT  70°F  FOR  NONLINEAR  MATERIAL  MODE! 


MATERIAL 

PROPERTY 

A 

B 

i] 

! 

o 

1 

~ Ert 

1.851 9x1 06psi 

.191829 

.314467 

Ezt 

1. 235x1 06psi 

.204557 

*293684 

vret 

.1 

.376247 

.169649 

vzet 

.1383 

.407264 

.0879388 

Grz 

,54xlOEpsi 

0 

1 

U„ 


1 psi 


Fig.  17,  Ref.  30,  ICO 
with  circumferential 
stress  gradient 
adjustment 

Fig.  14,  Ref.  30,  ICO 
Fig.  7,  Ref.  30,  16K9 
Fig.  14,  Ref.  30,  ICO 
Fig.  28,  Ref.  14,  Avg. 


tensile  and  compressive  material  properties  are  not  investigated  with  the  non- 
linear procedure.  Jortner's  across-grain  tension  results  in  Fig.  14  of  Ref. 

30  are  used  to  determine  the  material  property  constants  for  Ezt  and  vzet. 

His  with-grain  tension  results  in  Fig.  17  of  Ref.  30  are  used  to  establish 
the  constants  for  Er|.(=Eet).  Values  of  the  remaining  two  independent  material 
properties  are  determined  from  results  from  billets  other  than  1C0-15.  Mate- 
rial property  constants  for  vrQt  are  identical  to  those  in  the  16K9  data  base 
in  Table  8.  The  remaining  material  property,  Grzt,  is  assumed  to  have  a con- 
stant value  of  .54  x 10®  psi.  The  value  of  has  a negligible  effect  on 
the  predicted  strain  results  in  the  gage  sections  of  the  test  specimens. 

A correction  for  the  circunferential  stress  gradients  is  applied  to  the 
stress  values  In  Fig.  17  of  Ref.  30  before  the  material  property  constants  for 
Ert  are  evaluated.  Oorther  [9,  30]  suggests  that  the  circumferential  stress 
at  the  outer  surface  of  the  specimens  is  approximately  related  to  the  mean 
stress  used  as  the  ordinate  in  Fig.  17  of  Ref  30  by  the  expression 
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* , = Ko  (33) 

outer  surface  *Tnean 

where 

K = .96  + .01 

The  mean  (thin-wall)  circumferential  stresses  are  multiplied  by  .96  to  obtain 
an  estimate  of  the  circumferential  stresses  in  the  specimen  at  the  location 
where  the  strains  were  measured. 

The  new  nonlinear  material  model  is  a significant  improvement  over  the 
previously  used  material  model  for  the  two  biaxial  tension  stress  ratios  in- 
vestigated. The  percentage  deviation  between  test  data  and  predictions  for 
the  (1:1)  stress  ratio  is-29.4%  for  the  SAAS  II  material  model  and  +3.1%  for 
the  new  nonlinear  material  model.  Similar  correlation  is  found  for  the 
(1:1.26)  stress  ratio: .a  -32.0%  deviation  is  obtained  for  the  SAAS  II  non- 
linear model  and  a -2.6%  deviation  was  obtained  for  the  new  nonlinear  material 
model.  The  percentage  deviations  discussed  in  this  section  are  based  on  the 
experimental  point  closest  to  the  predicted  point.  For  both  biaxial  tension 
stress  ratios,  the  SAAS  II  material  model  predictions  are  a significant  under- 
estimate of  the  strain  response.  Based  on  the  strain  correlations,  the  new 
nonlinear  material  model  has  excellent  potential  for  description  of  the  bi- 
axial softening  phenomenon. 

Strain  predictions  with  the  new  nonlinear  material  model  for  the  (-.64:1) 
and  (-1:1)  stress  ratios  are  poorer  than  for  the  biaxial  tension  stress  ratios. 
However,  in  those  predictions,  the  multimodulus  character  of  the  material  is 
ignored!  The  deviations  between  the  strain  predictions  of  the  new  nonlinear 
material  model  and  the  closest  experimental  results  are  +15.7*  for  the  (-.64:1) 
stress  ratio  and  +18.0%  for  the  (-1:1)  stress  ratio,  (fcviously,  the  different 
moduli  of  ATJ-S  graphite  u.sder  tension  and  compression  must  be  accounted  for. 

The  apparently  good  agreement  (2.7%  deviation)  for  the  SAAS  II  material 
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model  is  both  coincidental  and  fortuitous  as  already  discussed.  That  is,  the 
SAAS  II  yield  criterion  is  obviously  crude  and  its  stress -strain  model  is  bi- 
linear instead  of  having  the  actual  nonlinear  behavior.  Those  facts,  coupled 
with  the  poor  predictions  for  the  SAAS  II  model  in  the  biaxial  tension  quad- 
rant, lead  to  a complete  invalidation  of  the  SAAS  II  (and  the  SAAS  III)  non- 
linear material  model  for  ATJ-S  behavior.  However,  the  new  nonlinear  material 
model  has  attractive  characteristics  in  biaxial  tension.  The  strain  correla- 
tions for  the  nonlinear  multiraodulus  procedures  are  examined  next. 


£.2.2  Strain  Correlations  for  Nonlinear  Multi modulus  MateriaLModel 


Predictions  with  the  16K9  Data  Base 


The  experimental  strains  obtained  at  roan  temperature  and  the  strains 
predicted  with  the  nonlinear  weighted  compliance  material  model  for  a maximum 
principal  thin-wall  stress  of  3550  psi  are  shown  in  Fig.  57.  lest  data  from 
four  different  billets  are  included  on  this  figure.  Because  a data  set  for 
each  billet  complete  enough  to  serve  as  a uniaxial  data  base  for  the.  nonlinear 
WCM  is  not  available,  the  16K9  data  base  in  Table  8 is  used.  Recall  that  this 
data  base  is  derived  from  a portion  of  the  results  of  Jortner's  uniaxial  off- 
-axis investigation  of  billet  15K9  described  in  Section  5.  Therefore,  the 
correlation  between  predicted  and  measured  strains  for  the  stress  ratios 
(1:1.26)  and  (0:1)  is  not  affected  by  billet-to-billet  variations  in  material 
properties. 

The  experimental  data  in  Fig.  57  are  from  Fig.  66  of  Ref.  8 end  Figs.  14- 
19  of  Ref.  30.  Because  of  friction  in  the  loading  apparatus  for  the  biaxial 
tubular  specimens,  the  measured  axial  load  is  in  error  by  less  than  1.51  [9, 
30].  Therefore,  the  axial  stress  in  the  specimen  is  approximately 
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FIGURE  57  PREDICTED  AND  MEASURED  STRAINS  IN  A HOLLOW  GRAPHITE  BIAXIAL 
TEST  SPECIMEN  - ROOM  TEMPERATURE  NONLINEAR  WCM  PREDICTIONS 
- 16K9  DATA  SET 


A similar  friction  effect  is  noted  for  the  compression  rod  results  in  the 
data  base  of  Table  8.  Thus,  a friction  correction  is  appropriate  for  the 
axial  stress  components  of  the  (1:0),  (1:1),  and  (1:1.26)  stress  states.  No 
friction  correction  need  be  applied  to  the  axial  compression  - internal 
pressure  biaxial  stresses  because  the  behavior  of  the  biaxial  specimens  is 
similar  to  that  of  the  uniaxial  compression  rods  from  which  the  stress-strain 
data  are  obtained.  That  is,  there  is  no  dissimilar  behavior  as  between  the 
uniaxial  tension  bar  and  the  biaxial  specimen  under  biaxial  tension. 

The  actual  experimental  strains  are  obtained  for  the  (1:0)  and  (1:1.26) 
stress  ratios  in  the  following  manner.  The  tensile  stress-strain  curve  in 
the  axial  direction  for  the  biaxial  specimen  is  entered  at  50  psi  less  than 
the  nominal  stress,  and  the  corresponding  strain  is  read.  (A  more  precise 
stress  correction  is  unrealistic  because  (1)  the  friction  error  is  approxi- 
mate and  (2)  50  psi  is  the  maximum  practical  resolution  on  the  plotted  stress- 
strain  curves.)  The  circumferential  strain  for  the  (1:1.26)  biaxial  stress  is 
read  directly,  (i.e.,  without  adjustment  for  friction).  The  two  strains  (axial 
and  circumferential)  are  then  plotted  in  Fig.  57  with  the  assumption  that  they 
are  for  the  same  biaxial  load.  The  experimental  strains  for  the  (1:1)  stress 
ratio  are  plotted  without  correction  because  stress-strain  curves  are  not  avail- 
able. Accordingly,  results  from  Fig.  66  of  Ref.  8 are  used  directly.  The  pre- 
dicted strains  are  obtained  by  use  of  the  nominal  stresses. 

The  predicted  strains  in  Fig.  57  are  obtained  with  the  MULTIW  computer  pro- 
gram described  in  Section  4.4.  For  nominal  stress  ratios  (1:1.26)  and  (-1:1), 
two  experimental  data  points  with  slightly  different  actual  stress  ratios  are 
presented.  Thus,  slightly  different  strain  predictions  are  obtained  for  each 
actual  stress  ratio.  Both  strain  predictions  for  the  (1:1.26)  stress  ratio 
nearly  coincide.  The  predictions  shown  for  the  (-1:1)  stress  ratio  are  com- 
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puted  for  the  actual  stress  ratio  of  the  billet  10V9  results  (inverted  tri- 
angle). Three  different  sets  of  strain  predictions  are  shown  for  each  of  the 
two  "mixed"  stress  ratios,  (-.64:1)  and  (-1:1).  The  different  strain  pre- 
dictions for  these  stress  ratios  are  a result  of  the  differences  in  the  three 
energy  approaches  investigated  (total  energy,  weighted  energy  with  n=l,  and 
divided  energy).  Strain  predictions  for  the  (-1:0)  stress  ratio  are  shown  in 
Fig.  57  despite  the  fact  that  no  strains  were  measured  on  biaxial  specimens 
at  that  stress  ratio  (the  biaxial  specimens  might  buckle).  These  predicted 
strains  for  the  hollow  graphite  specimen  are  based  on  uniaxial  stress-strain 
data  for  compression  rods.  The  strain  prediction  is  identical  for  all  energy 
approaches  with  the  nonlinear  RCM  for  the  (-1:0)  stress  ratio. 

Because  of  billet-to-billet  variation  in  material  properties,  the  cor- 
relations in  Fig.  57  are  mainly  qualitative  in  nature.  The  differences  be- 
tween predicted  and  measured  strain  values  for  the  (1:0)  stress  ratio  are  at- 
tributed mainly  to  differences  in  billet  stiffnesses.  The  predicted  strains 
are  based  on  the  16K9  billet  data  whereas  the  measured  strains  are  from  tests 
on  specimens  from  billets  ICO  and  10V9.  Based  on  the  strain  correlations  for 
the  (1:0)  stress  ratio  case,  billet  16K9  is  stiffer  in  the  across-grain  (z) 
direction  than  either  billets  ICO  or  10V9.  All  three  energy  approaches  in 
the  nonlinear  WCM  result  in  the  identical  strain  predictions  for  the  (1:0) 
loading  case.  Identical  strain  predictions  also  result  for  the  other  two 
uniaxial  loading  cases,  (0:1)  and  (-1:0),  and  the  two  biaxial  tension  loading 
cases,  (1:1)  and  (1:1.26). 

The  predicted  and  measured  strains  for  the  equal  biaxial  tension  stress 
ratio  case,  (1:1),  agree  to  within  5.5%.  Here,  again,  the  strain  correlations 
could  be  affected  by  billet-to-billet  variations  in  material  properties.  For 
the  two  (1:1.26)  cases,  the  deviations  between  predicted  and  measured  strains 
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are  4.6%  and  7.4%.  Recall  that,  for  these  cases  and  the  (0:1)  case,  the  data 
base  material  properties  and  the  measured  strains  are  obtained  from  tests  of 
specimens  from  the  same  billet.  Thus,  no  undesired  effects  are  introduced 
in  the  strain  correlations  by  variations  in  material  properties  between  bil- 
lets. However,  some  variations  between  predicted  and  measured  strain  results 
can  be  expected  due  to  the  gradient  in  circumferential  stress  through  the 
thickness  of  the  biaxial  test  specimens.  The  stresses  at  the  outer  surface 
of  the  specimens  where  the  strains  are  measured  can  be  different  from  the 
"thin-wan"  stresses  used  as  a basis  for  the  predictions.  The  circumferential 
stress  gradient  could  be  the  major  factor  in  the  7.7%  deviation  between  pre- 
dicted and  measured  strains  found  for  the  (0:1)  stress  ratio  case.  Consistent 
billets  are  not  used  to  obtain  the  predicted  and  measured  strains  for  the 
(-.64:1)  and  (-1:1)  stress  ratio  cases  in  Fig.  57.  Consequently,  an  evalua- 
tion of  the  predictions  for  the  different  energy  approaches  will  be  deferred 
until  Section  6.2.3. 

Strain  predictions  for  the  nonlinear  RCM  obtained  with  the  MULTI R com- 
puter program  and  strains  measured  in  room  temperature  tests  are  shown  in  Fig. 
58.  The  strain  predictions  are  for  a maximum  principal  thin-wall  stress  of 
3550  psi.  The  experimental  strains  presented  in  Fig.  58  are  identical  to  those 
shown  in  Fig.  57.  Material  property  set  Ezc  described  in  Section  5.3.1  was 
used  as  a uniaxial  data  base  for  the  nonlinear  restricted  compliance  material 
model.  The  material  property  constants  for  the  independent  material  proper- 
ties in  the  data  base  are  obtained  from  billet  16K9  values  in  Table  8.  Due 
to  the  composition  of  the  Ezc  property  set,  the  predictions  of  the  nonlinear 
RCM  in  Fig.  58  and  those  of  the  nonlinear  WCM  in  Fig.  57  are  identical  for 
the  biaxial  tension  and  uniaxial  tension  cases.  When  the  divided  energy  ap- 
proach is  used  in  the  nonlinear  RCM  with  the  Ezc  property  set,  the  transverse 
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strain  prediction  for  the  uniaxial  compression  stress  state  (-1:0)  is  differ- 
ent from  the  common  prediction  of  the  total  and  weighted  energy  approaches. 
With  the  divided  energy  approach,  the  implied  off-diagonal  compliances  are 
constant  on  successive  iterations  because  the  off-diagonal  compliances  are 
always  determined  from  the  linear  tensile  material  properties. 

The  strain  predictions  with  the  nonlinear  WCM  and  RCM  are  different  for 
the  mixed  biaxial  stress  states  in  Figs.  57  and  58.  The  per.entage  deviations 
in  the  strain  predictions  with  the  two  material  models  are  shown  in  Table  10. 


TABLE  10  PERCENTAGE  DEVIATIONS  OF  NONLINEAR 
RCM  STRAIN  PREDICTIONS  FROM  NONLINEAR  WCM  STRAIN  PREDICTIONS 


NOMINAL 

PERCENTAGE  DEVIATI 

DNS 

LOADING 

DIVIDED 

(V°c> 

ENERGY 

ENERGY 

ENERGY 

(-.64:1) 

2.4 

1.6 

1.5 

(-1:1) 

3.4 

2.3 

» 

CM 

The  strain  predictions  with  the  nonlinear  WCM  are  used  as  the  baseline  values 
for  calculation  of  percentage  deviations.  The  largest  percentage  deviation 
in  predicted  strains  is  3.4%.  Thus,  the  effects  on  strain  predictions  of  the 
differences  in  the  multi  modulus  formulations  of  the  two  material  models  are 
relatively  small  for  the  loading  coi.ditions  and  material  properties  investi- 
gated. The  agreement  between  predicted  and  measured  strains  found  for  one 
material  model  could  be  used  in  a general  assessment  of  both  material  models. 

6.2.3  Strain  Correlations  for  Nonlinear  Restricted  Compliance  Material  Model 
Predictions  with  the  ICO  Data  Base 

The  experimental  strains  and  the  strains  from  two  numerical  prediction 


115 


techniques  are  shown  in  Fig.  59.  Two  separate  sets  of  strain  predictions  for 
a maximum  principal  thin-wall  stress  of  3550  psi  are  shown  for  t1"  . onlinear 
restricted  compliance  material  model.  The  uniform  stress  model  in  the  MULTIR 
computer  program  is  one  prediction  technique.  The  355  element  model  of  the 
specimens  discussed  in  Section  6.2.1  is  used  with  a modified  version  of  the 
SAAS  III  computer  program.  Both  the  nonlinear  restricted  compliance  and  the 
nonlinear  weighted  compliance  material  models  with  all  three  energy  approaches 
are  incorporated  in  the  modified  version  of  the  SAAS  III  program.  Convergence 
of  the  nonlinear  material  models  with  changing  stresses  on  each  iteration  can 
be  studied  with  the  finite  element  procedure.  Moreover,  stress  and  strain 
gradients  can  be  predicted  in  the  specimens.  Strain  predictions  shown  for 
the  nonlinear  RCM  in  the  modified  SAAS  III  program  are  extrapolated  to  the 
outer  surface  of  the  specimens  where  the  experimental  strains  were  measured. 

The  experimental  results  shown  in  Fig.  59  for  all  but  one  case  are  ob- 
tained from  tests  of  specimens  from  a single  billet  (1C0-15)  in  contrast  to 
Figs.  57  and  58.  The  measured  strains  presented  for  the  equal  biaxial  ten- 
sion stre?  * state  were  found  in  tests  of  a specimen  from  billet  3R9-33. 

Tv/  reduce  effects  caused  by  differences  in  billet  stiffness,  specimen 
design,  and  test  procedures  in  the  evaluation  of  the  nonlinear  multimodulus 
material  models,  uniaxial  tests  of  the  tubular  specimens  were  used,  where 
possible,  to  establish  the  material  property  constants.  The  material  prop- 
erty constants  [A,  B,  C,  and  UQ  of  Eq.  (5)]  for  the  six  independent  material 
properties  in  the  nonlinear  RCM  property  set  Ezc  are  given  in  billet  coordi- 
nates in  Table  11.  These  material  property  constants  are  used  to  obtain  both 
sets  of  predictions  for  the  nonlinear  RCM.  The  origin  of  the  material  prop- 
erty constants  for  four  of  the  material  properties  in  this  table  (Ert,  Ezt, 
and  v26t)  which  are  repeated  from  Table  9 is  discussed  in  Section  6.2.1.  As 
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TABLE  11 

ICO  DATA  BASE  AT  70° F FOR  NONLINEAR  RCM 


MATERIAL 

PROPERTY 

A 

B 

Ert 

1.851 9x1 06psi 

.191829 

Ezt 

l.Z35xl 06psi 

.204557 

vret 

.1 

.376247 

vzet 

.1383 

.407264 

E;zt 

1. 385x1 06psi 

.132548 

m 

Nl 

o 

1.04xl06psi 

.114414 

DATA 

SOURCE 


Fig.  17,  Ref.  30,  ICO 
with  circumferential 
stress  gradient 
adjustment 

Fig.  14,  Ref.  30,  ICO 

Fig.  7,  Ref.  30,  16K9 

Fig.  14,  Ref.  30,  ICO 

Fig.  5,  Ref.  30,  16K9 

Fig.  9,  Ref.  30,  16K9 
Adjusted  for  billet 
ICO  stiffness 


was  the  case  with  v , material  property  constants  for  Erzt  are  obtained  from 
results  from  billet  16K9  instead  of  ICO.  Consequently,  the  material  property 
constants  for  vf0t  and  E'rzt  are  identical  in  both  the  16K9  (Table  8)  and  ICO 
(Table  11)  data  bases. 

Material  property  constants  for  Ezc  are  obtained  from  results  from  bil- 
let 16K9.  Uniaxial  stress-strain  results  are  not  available  for  the  tube 
specimens  from  billet  ICO  under  compression  loading  (the  specimens  would 
buckle! ).  Thus,  an  adjustment  for  the  apparent  difference  in  stiffnesses  of 
billets  16K9  and  ICO  is  made.  The  uniaxial  strain  results  in  Fig.  9 of  Ref. 

30  are  obtained  from  the  across-grain  compression  loading  of  the  rod  speci- 
mens from  billet  16K9.  At  a given  stress  level,  the  strain  from  Fig.  9 of 
Ref.  30  is  multiplied  by  the  ratio  of  the  strains  found  in  across-grain  uni- 
axial tension  tests  of  tube  specimens  from  billet  ICO  and  rod  specimens  from 
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billet  16K9.  The  resultant  strain  is  used  to  evaluate  the  material  property 
constants  for  Ezc. 

The  ICO  data  base  is  used  with  both  nonlinear  RCM  prediction  techniques 
considered  in  Fig.  59.  For  the  convergence  criterion  used  in  the  MULTIR  and 
modified  SAAS  III  programs  (relative  change  in  total  strain  energy  < .005), 
the  procedures  converge  in  four  to  six  iterations  for  each  of  the  cases  ana- 
lyzed. Recall  that  the  differences  in  predictions  of  the  two  nonlinear  RCM 
approaches  are  due  to  stress  gradients  in  the  specimens.  For  the  five  stress 
ratios  which  have  a thin-wall  circumferential  stress  of  3550  psi,  the  circum- 
ferential strains  predicted  with  the  355  element  grid  in  the  SAAS  III  program 
are  between  4%  and 6%  less  than  those  predicted  with  the  MULTIR  program  with 
the  same  energy  approach. 

Strains  predicted  with  the  nonlinear  RCH  agree  well  with  measured  strains 
for  the  two  biaxial  tension  stress  ratios  investigated  with  the  modified  SAAS 
III  approach.  For  the  equal  biaxial  tension  stress  ratio,  a 3%  deviation  be- 
tween predicted  and  measured  strains  is  found.  The  predictions  with  the  modi- 
fied SAAS  III  approach  for  the  nominal  (1:1.26)  stress  ratio  are  based  on  the 
stress  ratio  of  the  upper  of  the  two  data  points.  Those  predictions  are  7% 
less  than  the  experimental  strains.  If,  instead,  the  stress  ratio  of  the  lower 
experimental  point  were  used,  better  (but  not  significantly  better)  correlation 
would  be  obtained.  This  speculation  is  based  on  the  observation  that  the  two 
predictions  with  the  MULTIR  program  are  essentially  on  top  of  one  another;  a 
similar  result  is  expected  for  the  modified  SAAS  III  results. 

The  actual  ratio  of  axial  stress  to  circumferential  stress  for  the  two 
data  points  with  a nominal  (1:1.26)  stress  ratio  in  Fig.  59  is  slightly  smal- 
ler for  the  upper  point  than  for  the  lower  point.  When  the  circumferential 
stress  is  3550  psi  in  Fig.  16  of  Ref.  30,  then  the  axial  stress  in  Fig.  15  of 


Ref.  30  for  the  upper  experimental  data  point  in  Fig.  59  is  lower  than  the 
axial  stress  for  the  lower  point.  The  upper  point  has  the  largest  strains, 
but  corresponds  to  a lower  axial  stress  in  contradiction  to  the  expected 
stress-strain  behavior.  Thus,  the  differences  between  the  two  experimental 
points  for  the  nominal  (1:1.26)  stress  ratio  are  probably  an  indication  of 
the  data  scatter  and  not  a result  of  differences  in  actual  stress  levels. 

Given  the  scatter  expected  in  experimental  results,  a 7%  deviation  may  repre- 
sent good  agreement. 

For  the  (1:1.26)  stress  ratio,  the  strain  predictions  of  the  nonlinear 
material  model  shown  in  Fig.  56  are  not  identical  to  the  strain  predictions 
of  the  nonlinear  RCM  shown  in  Fig.  59  predominantly  because  of  the  different 
stress  ratios  used.  However,  even  for  the  identical  stress  ratios,  a small 
difference  between  the  two  sets  of  predictions  could  occur.  A radial  stress 
gradient  exists  through  the  thickness  of  the  specimen.  Thus,  a tri axial 
stress  state  must  actually  be  analyzed.  Different  material  properties  are 
used  to  determine  the  compliances  in  the  two  material  models  since  the  radial 
stress  is  a compressive  principal  stress. 

In  the  biaxial  loading  quadrant  with  mixed  tensile  and  compressive  load- 
ing, the  correlation  between  predicted  and  measured  strain  response  is  depen- 
dent on  the  energy  approach  used.  The  percentage  deviations  for  the  (-.64:1) 
data  point  and  one  of  the  (-1:1)  data  points  are  presented  in  Table  12.  These 
percentage  deviations  are  based  on  the  experimental  strain  results.  Results 
for  both  the  MULTIR  and  the  modified  SAAS  III  formulations  of  the  nonlinear 
RCM  are  presented.  Percentage  deviations  are  shown  for  the  MULTIR  program 
with  all  three  of  the  energy  approaches.  Only  the  weighted  energy  approach 
is  investigated  with  the  modified  SAAS  III  program.  The  actual  stress  ratios 
for  the  two  (-1:1)  data  points  are  (-.97:1)  and  (-.99:1).  The  percentage 
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TABLE  12 

PERCENTAGE  DEVIATIONS  OF  NONLINEAR  ROM 
STRAIN  PREDICTIONS  FROM  MEASURED  STRAINS 


deviations  in  Table  12  are  given  for  the  (-.99:1)  data  point.  To  determine 
the  effects  of  stress  gradients  in  the  specimens  on  the  predicted  strain 
response,  the  (-.99:1)  stress  ratio  is  investigated  with  the  nonlinear  RCM 
in  the  modified  SAAS  III  computer  program*  The  lower  of  the  two  (-1:1)  ex- 
perimental points  shown  in  Fig.  59  is  the  strain  result  for  the  (-.99:1) 
stress  ratio.  Thus,  the  strain  predictions  denoted  by  the  diamond  symbol  must 
be  compared  with  the  experimental  strains  for  the  lower  of  the  two  circles. 

Based  on  comparison  of  the  percentage  deviations  between  predicted  and 
experimental  strains  presented  in  Table  12,  the  weighted  or  divided  energy 
approaches  with  deviations  of  3-9%  are  superior  to  the  total  energy  approach 
with  deviations  of  20-24%.  A recommendation  as  to  which  approach,  the  di- 
vided energy  or  the  weighted  energy,  should  be  generally  used  cannot  be  made 
based  on  the  percentage  deviations  in  Table  12*  However,  from  the  standpoint 
of  the  theoretical  formulations  of  the  nonlinear  RCM,  the  weighted  energy  ap- 
proach might  be  preferable.  With  the  weighted  energy  approach,  all  material 
properties  are  evaluated  at  the  sane  energy  level.  In  that  manner,  the  addi- 
tional restrictive  implications  of  the  evaluation  of  tensile  and  compressive 
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material  properties  at  different  energy  levels  are  avoided. 


6.3  FRACTURE  STRAIN  CORRELATIONS 

The  experimental  stress-strain  results  obtained  by  Jortner*  [6-8J  in  the 
biaxial  fracture  investigation  at  room  temperature  are  used  to  further  evaluate 
the  nonlinear  multimodulus  material  models.  In  this  section,  measured  strains 
at  failure  are  compared  with  predicted  strains.  In  fact,  these  fracture  strain 
results  were  obtained  by  further  loading  of  the  same  specimens  discussed  in 
the  preceding  section.  Note  that  a failure  criterion  is  not  included  in  the 
modified  SAAS  III  computer  program  nor  in  MULTIW  or  MULTIR,  Thus,  only  the 
relationships  between  failure  stresses  and  strains  are  predicted  herein,  not 
the  stresses  and  strains  which  will  produce  failure.  The  failure  stress  lev- 
els are  input  quantities  used  to  predict  the  strains  at  failure. 

Hollow  test  specimens  (see  Fig.  54)  were  loaded  at  roan  temperature  to 
failure.  Stresses  and  strains  at  fracture  found  by  Jortner  [7]  for  six  of  the 
tubular  test  specimens  from  billet  16K9-27  are  presented  in  Table  13. 

The  predictions  with  both  the  nonlinear  WCM  and  the  nonlinear  RCM  are 
shown  along  with  measured  strains  in  Figs.  60  and  61,  respectively.  Data 
base  material  properties  in  Table  8 are  used  with  the  MULTIR  and  MULTIW  com- 
puter programs  to  obtain  the  predicted  strains  in  these  figures.  Recall  that 
the  material  property  constants  presented  in  Table  8 are  determined  from  uni- 
axial tests  of  specimens  from  billet  16K9-27.  Thus,  billet-to-billet  varia- 
tions in  material  properties  are  not  introduced  in  the  strain  predictions  in 
Figs.  60  and  61.  The  effects  of  stress  gradients  on  strains  in  the  hollow 
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specimens  are  expected  to  be  similar  in  nature  and  perhaps  slightly  larger 
than  those  found  in  the  strain  predictions  of  Section  6.2.3.  Moreover,  those 
stress  gradient  effects  can  be  evaluated  only  by  use  of  the  relatively  expen- 
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TABLE  13 

FRACTURE  STRESSES  AND  STRAINS  AT  70° f 
FROM  BILLET  16K9-27  {Ref.  7,  p.  52) 


live  finite  element  approach,  flag#  the  essential  characterise cs:  of  the 
stress  gradient  effects  are  known  well  enough  that  the  finite  element  analy- 
sis expense  is  not  justified.  Basically,  the  stress  gradient  is  believed  to 
be  the  main  cause  cf  differences  between  predicted  and  measured  strains  for 
the  specimen  under  internal  pressure. 

The  predictions  with  the  nonlinear  WCM  and  the  nonlinear  RCH  are  identi- 
cal for  all  stress  states  except  the  mixed  biaxial  stress  states  (b.  = -400 

^ ' ^|l  - 

psi , oQ  = 4890  psi ) for  the  loading  conditions  and  material  property  sets  in- 
vestigated. The  differences  in  predicted  strains  for  the  mixed  biaxial  stress 
states  are  small  and  not  even  perceptible  in  Figs.  60  and  61.  The  largest 
deviation  between  predictions  with  the  two  material  models  based  on  the  pre- 
diction of  the  nonlinear  RCH  is  approximately  .5%.  Approximate  percentage 
deviations  between  measured  strains  and  strains  predicted  with  the  nonlinear 
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TABLE  14 

PERCENTAGE  DEVIATIONS  OF  NONLINEAR  RCM 
STRAIN  PREDICTIONS  (WEIGHTED  ENERGY  APPROACH)  FROM  FRACTURE  STRAINS 


CASE 

NUMBER 


1 

2 

3 


°a 

(psi) 

°c 

(psi) 

PERCENTAGE 

DEVIATION 

3980 

0 

8.2 

3900 

0 

2.8 

3870 

4900 

10.4 

3550 

4490 

5.5 

0 

4710 

7.0 

-400 

4890 

16.2 

strains  are  used  as  base  values  in  the  computation  of  the  percentage  devia- 
tions. As  can  be  seen  in  Table  14,  the  smallest  deviation,  2.8%,  occurs  for 
the  uniaxial  stress  state  (oa  = 3900  psi,  oc  * 0)  and  the  largest  deviation, 
16.2%,  occurs  for  the  mixed  biaxial  stress  state.  The  percentage  deviations 
in  Table  14  for  cases  3-6  would  be  reduced  if  the  effects  of  stress  gradients 
were  evaluated.  The  overall  agreement  between  predicted  and  measured  strain 
results  is  only  fair  since  the  largest  deviation  is  16.2%.  However,  as  frac- 
ture is  approached,  the  stress-strain  curves  become  more  nonlinear  than  at 
3550  psi  maximum  principal  stress.  Thus,  fracture  strain  predictions  should 
be  expected  to  be  less  accurate  than  at  lower  stress  levels. 


6.4  2000 °F  STRAIN  CORRELATIONS 

The  2000°F  stress-strain  behavior  of  ATJ-S  graphite  under  multiaxial 
load1’  .gs  was  investigated  by  Jortner  [7-9].  His  measured  biaxial  strain  re- 
sults are  compared  with  the  predictions  of  the  nonlinear  RCM  in  this  section, 
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The  hollow  test  specimens  used  in  the  2000°F  biaxial  investigations  are 
identical  in  shape  to  those  used  in  the  room  temperature  investigations  (see 
Fig.  54).  A thin  nickel  bladder  was  used  inside  the  2000°F  specimens  to  con- 
tain the  internal  pressurizing  gas  instead  of  the  rubber  bladder  used  in  the 
room  temperature  specimens.  Strains  at  2000°F  were  measured  with  optical  ex- 
tensometers.  (Recall  that  strain  gages  were  used  in  room  temperature  tests.) 

Several  different  ratios' of  axial  stress  to  circumferential  stress  (oa:oc) 
were  investigated.  The  loading  directions  were  aligned  with  the  nominal 
principal  material  direction.  Details  of  the  specimen  design,  loading  appar- 
atus, and  test  procedures  are  found  in  Refs.  7 and  8. 

To  implement  the  nonlinear  multimodulus  strain  prediction  procedures, 
a uniaxial  data  base  is  required  to  establish  the  relations  between  material 
properties  and  total  strain  energy.  The  2000°F  experimental  biaxial  strain 

! 

results  were  obtained  by  Oortner  [8]  for  ATO-S  graphite  billet  1C0-15.  Un- 
fortunately, uniaxial  test  data  at  2000°F  are  not  available  for  billet  1C0-15. 

The  uniaxial  data  bases  used  for  the  strain  correlations  in  this  section  are 
obtained  from  the  average  results  of  Starrett  and  Pears  [14].  Those  results 
are  the  average  of  test  results  from  a large  number  of  tests  performed  on 
numerous  billets. 

Material  property  constants  for  the  "average"  2QQ0°F  data  base  are  pre- 
sented in  Table  15.  The  constants  for  the  material  properties  Ert,  Ezt,  and 
Ezc  are  derived  from  the  stress-strain  curves  in  Figs.  4,  8,  and  20  of  Ref. 

14,  respectively.  For  the  balance  of  the  material  properties  in  the  data 
set  in  Table  15,  stress-strain  data  are  not  available  to  define  material  prop- 
erty variation  with  total  strain  energy.  Thus,  constant  values  are  used  for 
vret»  vzet’  and  Erzt*  Some  raw  transverse  deformation  data  is  presented  by 
Starrett  [33],  However,  these  data  are  results  from  individual  tests  and  are 
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TABLE  15 

AVERAGE  DATA  BASE  AT  2000° F 


f 


MATERIAL 

~ 

~1 

DATA 

PROPERTY 

A 

B 

C 

uo 

SOURCE 

Ert 

1.88xl06psi 

.0321362 

. 788298 

1 psi 

Fig.  4,  Ref.  14 

Ezt 

1.37 5x1 06psi 

.0527850 

.665860 

Fig.  8,  Ref.  14 

vret 

.128 

0 

1 

Fig.  26,  Ref.  14 

vzet 

.108 

0 

1 

vrzt  and  comP^ance 

relation 

Erzt 

1.4875x10® 

0 

1 

Compliance  transfor- 
mation relationship 
with  Grz  from  Fig. 
28,  Ref.  14 

Ezc 

1.58xl06psi 

.156581 

.352135 

Fig.  20,  Ref.  14 

t 

-1 

1 


not  compatible  with  “average"  results.  The  material  property  constant  A for 
the  direct  modulus  at  45°  to  the  r-z  principal  material  directions  is  esti- 
mated from  the  value  of  Grz  in  Fig.  28  of  Ref.  14  and  a compliance  transforma- 
tion relationship. 

A data  base  suitable  for  the  nonlinear  WCM  was  not  developed  because  the 
uniaxial  stress-strain  data  are  inadequate.  Strains  are  correlated  only  with 
the  nonlinear  RCM.  Differences  between  the  predictions  of  the  nonlinear  RCM 
and  the  nonlinear  WCM  are  small  in  the  room  temperature  investigation  in  Sec- 
tion 6,2.2.  In  fact,  for  the  material  property  sets  used  in  Section  6.2.2, 
the  predictions  with  the  nonlinear  RCM  and  the  nonlinear  WCM  are  identical 
for  the  biaxial  tension  stress  states.  An  RCM  property  set  with  the  same 
composition  of  independent  material  properties  is  used  in  this  section.  Ac- 
cordingly, the  results  for  the  nonlinear  RCM  might  be  representative  of  those 
expected  for  the  nonlinear  WCM. 

The  experimental  strains  obtained  at  2000°F  and  the  strains  predicted 
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with  the  nonlinear  restricted  compliance  material  model  are  shown  in  Fig.  62. 
The  experimental  data  points  in  Fig.  62  for  a maximum  principal  stress  of 
3550  psi  are  obtained  from  Fig.  67  of  Ref.  8.  The  average  data  base  in  Table 
15  is  used  with  the  MULTIR  computer  program  to  obtain  the  predicted  values. 

The  strain  correlations  for  the  uniaxial  loading  cases  (1:0)  and  (0:1)  are 
evidence  that  the  1C0-15  billet  stiffnesses  differ  from  the  average  values. 

In  the  axial  (across-grain)  direction,  billet  1C0-15  is  probably  not  as  stiff 
as  the  "average"  billet.  The  differences  between  predicted  and  measured 
strains  for  the  (0:1)  case  is  possibly  due  in  part  to  stress  gradients  in  the 
tube  specimens. 

To  obtain  a more  effective  evaluation  of  the  three  biaxial  loading  con- 
ditions, the  predicted  strains  for  the  two  uniaxial  loading  conditions  in 
Fig.  62  are  forced  to  approximately  correspond  to  experimental  results  by  ad- 
justments of  the  uniaxial  data  base  material  properties.  The  values  of  mate- 
rial property  constant  A for  three  of  the  six  independent  material  properties 
are  altered  in  the  following  manner.  The  values  of  A for  Ert  and  Ezt  are 
scaled  by  the  ratios  of  the  apparent  experimental  stiffness  at  3550  psi  to 
the  stiffness  predicted  with  the  average  data  base.  An  average  of  experimen- 
tal strains  is  used  to  define  the  apparent  experimental  stiffness.  For 
Ert  (=EQt ) , a correction  for  stress  gradient  is  also  made.  The  thin-wall 
stress  is  multiplied  by  .97  to  obtain  an  estimate  of  the  stress  at  the  outer 
surface  of  the  specimen  where  the  strains  were  measured.  The  factor  .97  is 
based  on  the  stress  gradient  predicted  for  the  tube  specimen  under  the  same 
uniaxial  loading  at  room  temperature.  The  new  value  for  vZQt  is  computed  from 
the  average  experimental  transverse  strain  for  the  (0:1)  stress  ratio  case. 

The  revised  2000°F  data  base  is  presented  in  Table  16.  Good  correlation  is 
ensured  at  a specific  strain  energy  level  by  these  adjustments.  However, 
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TABLE  16 

REVISED  DATA  BASE  AT  2000° F 


MATERIAL 

PROPERTY 

A 

B 

C 

H 

DATA 

SOURCE 

Ert 

2.01 35x1 06ps i 

.0321362 

. 788298 

1 P 

si 

Revised  Fig.  4,  Ref.  14 

E2t 

1.2495xl06psi 

.0527850 

.665860 

Revised  Fig.  8,  Ref.  14 

vr0t 

.128 

0 

1 

Fig.  26,  Ref.  14 

vzet 

.093 

0 

1 

Estimated  from  uniaxial 
loadings  herein  on  Fig. 
59 

Erzt 

1. 4875x1 06psi 

0 

1 

Compliance  transforma- 
tion relationship  with 
Grz  from  Fig.  28,  Ref.  14 

Ezc 

1.4357xl06psi 

.156581 

.352135 

< 

Revised  Fig.  20,  Ref.  14 

strain  correlation  at  other  energy  levels  is  still  an  unresolved  question.  In 
any  event,  the  available  experimental  results  are  not  adequate  for  more  re- 
fined adjustments. 

Strain  predictions  with  the  nonlinear  RCM  for  the  revised  2000° F data 
base  are  presented  in  Fig.  63  along  with  experimental  strains  at  a maximum 
principal  thin-wall  stress  of  3550  psi.  The  predicted  strains  for  the  (0:1) 
stress  ratio  case  do  not  coincide  with  the  average  experimental  strains  be- 
cause the  effects  of  stress  gradients  cannot  be  predicted  with  the  MULTIR  com- 
puter program  used  to  obtain  all  the  predicted  strain  results.  For  the  two 
biaxial  tension  stress  ratios,  (1:1.10)  and  (1:1.35),  the  predicted  circum- 
ferential strains  are  too  small  possibly  because  of  the  assumption  of  a con- 
stant tensile  Poisson's  ratio  (vz6t)  independent  of  the  strain  energy  level. 
However,  if  v2Qt  did  decrease  with  increasing  energy  level,  then  the  axial 
strains  as  well  as  the  circumferential  strains  would  increase. 
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FICUBE  63  PREDICTED  AT©  MEASURED  STRAINS  IN  A HOLLOW  GRAPHITE  BIAXIAL 
TEST  SPECIMEN  - NONLINEAR  RCM  PREDICTIONS  - 2000#F  REVISED 
DATA  SET 


A quantitative  discussion  of  the  strain  correlation  might  be  inappro- 
priate for  the  2000° F results  because  of  possible  inaccuracies  in  the  experi- 
mental results.  Jortner  [9]  states  that  the  use  of  the  optical  extensometers 
might  induce  "relatively  large  experimental  uncertainties"  in  the  measured 
strain  results.  He  also  states  that  the  nickel  bladder  is  so  stiff  that  it 
cannot  be  ignored  in  comparison  to  the  specimen  stiffness  (as  is  done  in  the 
present  correlations).  Even  with  the  uncertainties  in  both  the  predicted  and 
measured  strains,  the  general  trend  of  the  predictions  does  agree  with  the 
2000° F experimental  results. 
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6.5  SUMMARY 

Biaxial  strain  predictions  with  the  various  nonlinear  material  models 
described  in  this  report  are  compared  with  several  types  of  experimental  data 
in  this  section.  The  comparisons  with  room  temperature  data  are  for  biaxial 
tension  and  mixed  tension  and  compression  at  a maximum  principal  thin-wall 
stress  of  3550  psi  as  well  as  at  higher  stresses  at  fracture.  Comparisons  are 
| also  made  at  elevated  temperature  (2000°F)  to  round  out  the  comparisons  of 

model  predictions  with  representative  nosetip  conditions. 

The  room  temperature  biaxial  tension  predictions  are  within  about  3%  of 
the  experimental  results.  Thus,  the  basic  relationship  between  the  secant 
moduli  of  the  nonlinear  stress-strain  curves  and  the  strain  energy  is  shown  to 
be  accurate.  Accordingly,  the  phenomenon  of  softening  under  biaxial  tension 
is  now  predictable  with  the  new  material  model. 

The  extensi«..is  of  the  new  material  model  to  different  room  temperature 
behavior  ir.  tension  and  expression  are  not  quite  as  accurate  as  in  biaxial 
tension.  The  strain  predictions  are  fairly  good  (within  9%),  however,  when 
the  weighted  energy  or  divided  energy  approaches  are  used  in  the  nonlinear  RCM 


as  opposed  to  the  total  energy  approach.  The  predictions  with  the  nonlinear 
WCM  are  slightly  different  from  those  with  the  nonlinear  RCM  for  the  loading 
and  material  properties  considered.  However,  uniaxial  stress-strain  data  are 
often  not  sufficient  to  apply  the  nonlinear  WCM. 

The  room  temperature  predictions  of  strains  at  fracture  are  less  accurate 
than  the  strain  predictions  at  3550  psi  maximum  principal  stress.  As  fracture 
is  approached,  the  stress-strain  curves  become  more  nonlinear  than  at  3550  psi. 
Thus,  the  larger  strain  prediction  inaccuracy,  up  to  10%  in  biaxial  tension 
and  16%  in  mixed  tension  and  compression,  should  be  expected. 

The  2000°F  strain  predictions  are  the  least  accurate  of  all  comparisons 
made.  However,  the  2000°F  comparisons  are  adversely  influenced  by  the  fact 
that  the  stiffness  of  the  pressure  retaining  bladder  inside  the  graphite  spec- 
imen is  not  accounted  for  in  the  present  analysis. 
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1 7.  CONCLUDING  REMARKS 
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New  material  models  for  deformation  behavior  of  artificial  graphite  under 
initial  loading  are  discussed.  These  new  material  models  are  based  on  a new 
deformation  theory  of  orthotropic  plasticity  and  can  be  used  in  finite  element 
computer  programs.  One  material  model  is  used  to  predict  nonlinear  stress- 
strain  behavior  of  ATO-S  graphite  under  arbitrary  uniaxial  and  biaxial  tension 
stress  states.  Two  other  material  models  are  investigated  for  prediction  of 
both  nonlinear  and  multimodulus  deformation  behavior.  The  same  nonlinear 
approach  is  common  to  all  three  new  material  models.  The  distinction  between 
the  two  nonlinear  multi modulus  material  models  is  based  on  the  multimodulus 
approach  used  or,  more  specifically,  on  the  determination  of  the  compliance 
matrix. 

Each  independent  material  property  in  the  nonlinear  approach  can  vary  in 
an  arbitrary  manner  as  actually  occurs  in  uniaxial  tests  and  is  related  to  the 
multiaxial  state  of  stress  and  strain  by  a strain  energy  function.  In  linear 
multimodulus  material  models,  the  composition  of  the  compliance  matrix  is  de- 
termined from  the  signs  of  the  principal  stresses.  In  the  nonlinear  multimod- 
ulus material  model,  the  selection  of  material  properties  and,  consequently, 
the  stress-strain  relations,  is  based  on  both  the  signs  of  the  principal 
stresses  and  the  magnitude  of  an  energy  function.  An  iteration  procedure  is 
used  in  the  nonlinear  material  model  to  simultaneously  satisfy,  to  a specified 
tolerance,  the  nonlinear  stress-strain  relations  and  the  material  property 
versus  strain  energy  equations.  An  extension  of  the  iteration  approach  is 
used  in  the  nonlinear  multimodulus  material  models  to  simultaneously  satisfy 
the  constraints  of  both  the  nonlinear  and  multimodulus  problems.  The  iteration 
procedures  are  stable  and  converge  in  six  iterations  or  less  for  a majority 
of  the  cases  examined. 
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The  practical  use  of  the  new  material  models  is  no  more  difficult  than 
the  material  models  of  Jones  and  Crose  [24,  32]  which  are  now  used  in  reentry 
vehicle  nosetip  stress  analysis.  For  their  models,  constants  are  determined 
for  the  best-fit  bilinear  representation  of  uniaxial  stress-strain  data.  In 
both  the  new  nonlinear  material  model  and  the  new  nonlinear  multimodulus  mate- 
rial models,  constants  are  required  for  the  equations  relating  material  prop- 
erties to  the  strain  energy.  The  necessary  material  property  representation 
can  be  made  by  the  people  who  measure  the  material  properties.  For  example, 
Starrett  and  Pears  [14]  display  both  the  experimental  stress-strain  curves  and 
the  bilinear  approximation.  They  could,  with  little  extra  effort,  present  simi- 
lar recommended  material  property  values  for  the  new  material  models.  The 
number  of  independent  material  properties  required  to  apply  the  nonlinear  mul- 
timodulus  material  models  is  dependent  on  the  multimodulus  procedure  used  as 
well  as  on  the  symmetry  of  the  body  and  loading  and  the  degree  of  material 
anisotropy.  In  the  restricted  compliance  material  model,  restrictions  between 
the  material  properties  determined  in  uniaxial  tension  and  uniaxial  compression 
are  imposed  to  maintain  compliance  matrix  symmetry  under  transformation  of  co- 
ordinates. As  a consequence,  the  number  of  independent  material  properties  is 
reduced.  Because  of  these  relations,  restrictions  are  placed  on  the  stress- 
strain  behavior  that  can  be  treated  with  the  restricted  compliance  approach. 

Not  all  materials  have  properties  for  which  the  restrictions  are  satisfied. 

In  particular,  the  Poisson's  ratios  may  decrease  under  increasing  tension  load, 
but  increase  under  increasing  compression  load.  Thus,  the  properties  for  each 
material,  to  which  the  model  is  to  be  applied,  must  be  examined  to  see  if  the 
restricted  compliance  model  is  applicable.  If  those  restrictions  are  not  met, 
then  another  model  must  be  used. 

The  application  and  evaluation  of  the  new  material  modejs  is  hampered  by 
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the  lack  of  suitable  uni  axial  stress-strain  data  to  determine  material  prop- 
erty constants  in  the  general  nonlinear  stress -strain  relations.  Poisson's 
ratio  data  are  inadequate  to  properly  establish  the  variation  with  strain 
energy.  All  complexities  of  data  acquisition  are  multiplied  by  two  when  mul- 
timodulus  materials  are  considered.  That  is,  all  data  must  be  gathered  wider 
tension  as  well  as  under  compression  loading.  Moreover,  a new  parameter,  the 
Young's  modulus  at  45®  to  principal  material  directions,  must  be  measured  in 
both  tension  and  compression.  These,  new  moduli  are  the  replacement  for  the 
usual  shear  modulus. 

The  validation  of  the  new  material  models  would  be  hastened  if  uniaxial 
and  multiaxial  test  data  were  measured  on  specimens  from  the  same  billet. 
Accordingly,  billet-to-billet  variations  in  material  properties  could  be 
eliminated  as  one  cause  for  disagreement  between  predicted  and  measured 
strains.  Then,  the  predictions  and  measurements  could  be  compared  without 
undesirable  influences  and  hence  in  the  most  scientific  manner. 

The  material  property  data  and,  in  particular,  the  material  property 
constants  in  the  new  material  models  are  reported  to  more  (usually  six)  sig- 
nificant figures  than  necessary.  The  accuracy  of  the  model  will  not  suffer 
ac  all  if  the  number  of  significant  figures  is  reduced  to  three.  In  fact, 
use  of  more  than  three  significant  figures  in  reentry  vehicle  nosetip  stress 
analysis  problems  is  not  justified  when  the  inescapable  variability  of  graph- 
ite is  considered. 

In  the  nonlinear  single  modulus  and  the  nonlinear  muitimodulus  material 
models,  several  energy  functions  are  investigated  for  relating  multiaxial  to 
wiiaxial  deformation  behavior.  Total  strain  energy  and  distortional  strain 
energy  are  investigated  for  the  nonlinear  single  modulus  material  model.  How- 
ever, application  of  distortional  strain  energy  to  anisotropic  materials  is 
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hindered  by  the  interdependence  of  the  distortional  and  dilatational  energy 
components.  Thus,  the  usual  distortional  energy-based  theories  for  isotropic 
metals  are  not  applicable  to  anisotropic  materials  such  as  ATJ-S  graphite. 

For  the  nonlinear  multimodulus  models,  three  energy  functions  are  studied: 

(1)  total  strain  energy,  (2)  divided  strain  energy  (the  strain  energy  of  ten- 
sile stresses  is  separated  from  that  due  to  compressive  stresses),  and  (3) 
weighted  strain  energy  (the  strain  energy  used  in  the  material  property-energy 
relations  is  a weighted  combination  of  the  tensile  and  compressive  strain  ener- 
gies). The  tensile  and  compressive  strain  energies  are  separated  to  describe 
the  different  deformation  behavior  under  tensile  and  compressive  loads.  How- 
ever, in  the  divide^  energy  approach,  the  off-diagonal  compliances  must  be 
constant  for  all  energy  levels  so  even  more  material  behavior  restrictions  are 
imposed  than  with  the  basic  restricted  compliance  model. 

Two  nonlinear  multimodulus  material  models  are  used  in  prediction  of  uni- 
axial deformation  behavior  of  ATJ-S  graphite  for  loads  applied  at  various  an- 
gles from  principal  material  directions  (off-axis  tests).  When  these  predic- 
tions are  compared  with  the  strain  response  measured  by  Jortner  [30],  the  best 
overall  • irrelation  exists  for  the  nonlinear  weighted  compliance  model.  In 
that  model,  the  tension  and  compression  compliances  are  weighted  in  accordance 
with  the  proportion  of  tensile  and  compressive  stresses  to  obtain  a symmetric 
compliance  matrix.  Differences  between  measured  strains  and  the  predictions 
with  the  nonlinear  restricted  compliance  material  model  are  attributed  mainly 
to  the  restrictions  on  tension  and  compression  compliances  inherent  to  the 
model.  That  is,  the  stress-strain  behavior  of  ATJ-S  graphite  does  not  closely 
fit  the  restrictions  on  tension  and  compression  compliances.  Different  strain 
predictions  are  obtained  when  mainly  tension  compliances  are  used  and  compres- 
sion compliances  are  generated  by  use  of  the  restrictions  than  when  mainly 
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compression  compliances  are  used.  At  5000  psi,  deviations  of  8 to  12%  in 
axial  strain  pt  - cions  occur  between  the  two  types  of  predictions.  Accord- 
ingly, the  nonlinear  weighted  compliance  material  model  is  preferable  to  the 
nonlinear  restricted  compliance  material  model  for  graphite. 

Excellent  correlation  is  obtained  between  predicted  and  measured  strains. 
With  the  exception  of  two  low  stress  level  predictions,  differences  between 
predicted  and  measured  strains  are  less  than  8%.  Thus,  use  of  the  general 
material  property  equation,  Eq.  (5),  and  the  simultaneous  solution  procedure 
for  material  property  constants  B and  C results  in  an  effective  description 
of  the  materia!  property-strain  energy  variations  for  both  tensile  and  compres- 
sive properties  of  ATJ-S  graphite  at  70°F.  The  form  of  the  present  equation 
may  be  generally  applicable  to  a wide  variety  of  engineering  materials.  Other 
forms  for  the  general  material  property-strain  energy  equation  could  be  inves- 
tigated. Alternative,  but  more  complicated,  expressions  might  result  in  even 
batter  representation  of  material  property  variations  with  strain  energy. 

Biaxial  strain  predictions  are  compared  with  strain  values  from  Jortner‘s 
tests  of  ATJ-S  graphite  [7-9,  30]  to  evaluate  the  new  material  models.  The 
new  nonlinear  material  model  is  a significant  improvement  over  the  previously 
used  material  model  for  the  two  biaxial  tension  stress  ratios  investigated 
with  both  material  models.  Strain  predictions  with  the  nonlinear  material 
model  due  to  Jones  and  Crose  [24,  32]  are  approximately  30%  lower  than  mea- 
sured strains.  For  the  new  nonlinear  material  model,  the  theoretical  results 
are  within  about  3%  of  the  measured  strains  in  biaxial  tension.  Thus,  the 
new  nonlinear  material  model  is  appropriate  for  description  of  the  biaxial 
softening  phenomenon. 


The  strain  correlations  for  the  nonlinear  restricted  compliance  model 
with  either  the  weighted  energy  or  divided  energy  approach  are  within  9%  for 
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biaxial  loading  cases  with  mixed  tensile  and  compressive  principal  stresses  at 
room  temperature.  Strains  are  overpredicted  by  about  20%  with  the  total 
strain  energy  approach.  Based  on  the  results  of  some  preliminary  strain  cor- 
relations, the  predictions  with  the  two  nonlinear  multimodulus  material  models 
are  only  slightly,  if  at  all,  different  for  the  biaxial  stress  states  investi- 
gated. Thus,  similar  strain  correlations  are  anticipated  for  the  nonlinear 
weighted  compliance  model. 

When  measured  strains  at  fracture  are  compared  with  strain  predictions 
from  the  nonlinear  multimodulus  material  models,  the  overall  agreement  is  only 
fair  (the  largest  deviation  is  approximately  16%).  The  measured  failure 
stress  levels  are  used  to  predict  the  strains  at  failure.  Measured  biaxial 
strain  results  for  tests  at  2000°F  are  also  compared  with  predictions  from 
the  nonlinear  restricted  compliance  model.  Uncertainties  exist  in  the  pre- 
dictions due  to  neglect  of  the  stiffness  of  a reinforcing  bladder  inside  the 
graphite  specimens.  Moreover,  uncertainties  also  exist  in  the  measurement 
of  strains  with  optical  extensometers . Even  with  uncertainties  in  both  the 
predicted  and  measured  strains,  the  general  trend  of  the  predictions  does 
agree  with  the  2000°F  experimental  results. 

The  experimental  biaxial  results  are  few  in  number  and  limited  to  only 
two  quadrants  of  biaxial  loading.  Thus,  a comprehensive  validation  of  the 
nonlinear  multimodulus  material  models  must  wait  future  experimental  investi- 
gation. Study  of  ATJ-S  graphite  behavior  in  the  biaxial  loading  quadrants 
with  (1)  axial  and  circumferential  compression  and  (2)  axial  tension  a.id  cir- 
cumferential compression  is  recommended.  Without  such  a study,  the  possible 
hardening  of  graphite  under  biaxial  compression  (as  opposed  to  softening  under 
biaxial  tension)  cannot  be  confidently  predicted.  Both  stress  states  exist 
simultaneously  in  reentry  vehicle  nosetips  s?  they  must  be  well  understood. 
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The  nonlinear  weighted  complianct  material  model  with  the  weignted  energy 
approach  is  recommended  for  use  in  graphite  reentry  vehicle  nosetip  stress  anal- 
ysis. This  recommendation  is  based  on  the  fact  that  better  strain  correlation 
exists  for  that  model  than  for  the  nonlinear  restricted  compliance  material 
model.  The  lack  of  success  with  the  latter  model  is  probably  due  to  the  fact 
that  the  compliance  restrictions  are  not  satisfied  for  graphite  material  prop- 
erties. The  nonlinear  weighted  compliance  material  model  should  be  regarded 
as  an  engineering  approach  to  the  nonlinear  deformation  problem.  The  various 
energy  approaches  (weighted  energy,  etc.}  fa!1  in  the  same  category  as  opposed 
to  the  more  scientific,  but  unsubstantiated,  approach  of  the  restricted  compli- 
ance material  model. 

Application  of  the  new  material  models  to  stress  analysis  of  engineering 
materials  other  than  artificial  graphites  should  be  investigated.  Further 
theoretical  refinements  could  be  considered  to  increase  the  applicability  of 
the  new  material  models.  The  energy  function  for  the  multiaxial  stress  states 
is  one  phase  of  the  material  models  which  might  be  improved.  The  relations 
between  tension  and  compression  compliances  in  the  nonlinear  restricted  com- 
pliance material  model  could  be  extended  to  arbitrary  orientations  of  the 
principal  material  directions  in  the  body.  Currently,  one  principal  mate- 
rial direction  must  be  aligned  with  a body  coordinate  direction.  Extension 
of  the  material  models  to  analysis  of  general  bodies  under  arbitrary  loads  is 
recommended.  At  present,  the  material  models  are  limited  to  analysis  of  axi- 
symmetrically  loaded  axisymmetric  bodies  composed  of  orthotropic,  transversely 
orthotropic,  or  isotropic  materials. 
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APPENDIX  A COMPARISON  OF  THE  TOTAL  STRAIN  ENERGY  APPROACH 
WITH  THE  DISTQRTIONAL  ENERGY  APPROACH  FOR  ISOTROPIC  MATERIALS 
The  distortional  energy  approach  is  used  to  determine  stresses*  strains, 
and  displacements  in  isotropic  materials  which  deform  nonlinearly  [34].  The 
distortional  energy  is  the  single  parameter  used  to  express  the  level  of  deform- 
ation attained.  Alternatively,  the  distortional  energy  can  be  related  to  a 

stress  intensity,  o-,  where 

1 1/2 

°i  = K2  + °y2  ' °x°y  + 3Txy}  (35) 

The  constant  of  proportionality  between  the  distortional  energy  and  the  stress 

intensity  is  1/(66)  or  (l+v)/(3E).  This  theory  is  also  called  deformation 

theory  of  plasticity. 

In  J2  theory,  an  equal  value  of  is  achieved  for  the  stress  state  ox=o, 
0^=0,  ana  txy=0  as  for  cx=Q,  Oy=0,  and  -rxy=  o//~3.  For  example,  stress  states 
with  cx=  72  ksi  and  Txy=  41.57  ksi  have  the  same  distortional  energy  and,  hence, 
have  directly  relateable  moduli  Esec  and  Gsec.  That  is,  a relation  between  Esec, 

Jsec’ 

( 36) 


v , and  Gcor , 


E 

6 = sec._ 

sec  2(1  ♦ v) 


is  enforced  at  all  time  where  v varies  with  the  deformation  level  or  energy  level. 
The  distortional  energy  is  used  first  to  calculate  Esec»  ^en Esec  used  t0  cal" 
culate  v,  and  finally  Gsec  is  calculated  from  Eq.  (35). 

However,  in  the  present  total  strain  energy  approach,  the  total  strain 
energy,  not  just  the  distortional  energy  component,  is  used.  First,  observe  that 
the  distortional  energy  is  only  a component  of  the  total  strain  energy  in  the  uni- 
axial stress-strain  curve  of  Fig. 64  (a).  On  the  other  hand,  the  distortional 
energy  is  equal  to  the  total  strain  energy  in  the  shear  stress-shear  strain  curve 
of  Fig. 64  (b).  In  the  total  strain  energy  approach,  Esec  is  predicted  first  even 


if  only  is  applied  (and  hence  only  distortional  energy  exists  for  isotropic 
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FIGURE  64  DISTORTIONAL  VERSUS  TOTAL  STRAIN  ENERGY 
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materials).  Given  a fixed  applied  stress  t,  we  use  the  elastic  moduli  to  find 
an  initial  strain  approximation.  The  associated  energy  (which  is  pure  distor- 
tional  energy)  is  used  in  the  relation  for  E$ec 


Esec  ■ A[1'B(IJ„)C] 


(37) 


to  obtain  a value  of  Esec  which  is  less  than  the  initial  value  A.  The  coeffi- 
cients A,  B,  and  C are  determined  by  fitting  a curve  to  a plot  of  E$ec  versus 
U for  the  uniaxial  stress  case  oy  versus  ex.  Thus,  for  a given  value  of  U,  a 
different  correspondence  exists  between  ox  and  tXy  than  in  J2  distortional  energy 
theory.  We  will  try  to  trace  this  correspondence  through  several  iterations  to 
predict  the  type  of  results  we  will  obtain  with  our  theory  in  comparison  to  02 
theory. 

First,  we  use  elastic  properties,  namely  G,  to  obtain  a value  of  U which 
happens  to  be  distortion  in  Fig.  65  . That  value  of  U is  then  used  in  the  re- 
lation for  Esec,  Eq.  (37).  However,  that  value  of  U is  obviously  greater  than 
distortion  for  the  uni  axial  case  ox=o  in  Fig. 65  indicated  with  point  A.  The 
corresponding  point  B for  J2  theory  is  shown  in  Fig. 65  . The  horizontally  cross- 
hatched  area  associated  with  point  B is  the  ^distortion  of  the  original  t-y 
curve  and  the  ^ of  the  o-c  curve  to  point  A.  Thus,  at  every  iteration,  the 
value  of  E$ec  is  consistently  higher  for  our  total  strain  energy  than  for  J2 
distortional  energy  theory.  Accordingly,  the  value  of  Gsec  is  higher  for  our 
theory  than  for  J2  theory.  Then,  because  Y=t/Gsec,  the  predicted  shear  strains 
for  our  theory  are  smaller  than  for  J2  theory  when  t is  specified.  Finally,  the 
total  strain  energy  and  simultaneously  the  distortional  strain  energy  predicted 
is  smaller  in  our  theory  than  in  J2  theory. 

Suppose  we  attempt  to  reverse  the  sequence  of  events  and  start  with  a pre- 
diction of  G$ec  instead  of  Esec.  First,  we  must  find  v in  terms  of  Gsec.  We 
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know  that 


• • i - - vf* 


but  because  cf  Eq.  (36), 

1 ,,  2(1+v)Gsec 

2 7 2(l+ve)fi 

wnereupon,  if  we  solve  for  v, 

I - K G,  /G 
_ 2 sec 

v - r 

1 + K G /G 
sec 


where 


= 2 ~ve 
1 


Equation  (40)  reduces  to  ve  if  Gsec=G.  Then,  knowing  Gsec  and  u,  we  can  predict 
Esec‘  The  ener9-v  used  t0  define  the  shear  stress-shear  strain  curve  is  only  dis- 
tortional  energy.  Therefore,  when  the  total  strain  energy  is  found  for  an  arbi- 
trary stress  state  (for  example,  the  uniaxial  case,  ox  versus  ex)  and  then  used 
to  predict  Gsecs  a larger  value  of  energy  is  used  than  in  distortional  energy 
theory.  Accordingly,  G^  is  smaller  for  this  total  strain  energy  approach  than 
in  distortional  energy  theory. 

All  of  the  foregoing  rationale  has  been  substantiated  by  numerical  results. 
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